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Abstract 

We study higher-dimensional interlacing Fibonacci sequences, generated via 
both Chebyshev type functions and m-dimensional recurrence relations. For 
each integer to, there exist both rational and integer versions of these sequences, 
where the underlying prime congruence structures of the rational sequence de¬ 
nominators enables the integer sequence to be recovered. 

From either the rational or the integer sequences we construct sequences 
of vectors in Q m , which converge to irrational algebraic points in R m . The 
rational sequence terms can be expressed as simple recurrences, trigonometric 
sums, binomial polynomials, sums of squares, and as sums over ratios of powers 
of the signed diagonals of the regular unit n-gon. These sequences also exhibit 
a “rainbow type” quality, and correspond to the Fleck numbers at negative 
indices, leading to some combinatorial identities involving binomial coefficients. 

It is shown that the families of orthogonal generating polynomials defining 
the recurrence relations employed, are divisible by the minimal polynomials 
of certain algebraic numbers, and the three-term recurrences and differential 
equations for these polynomials are derived. Further results relating to the 
Christoffel-Darboux formula, Rodrigues’ formula and raising and lowering op¬ 
erators are also discussed. Moreover, it is shown that the Mellin transforms of 
these polynomials satisfy a functional equation of the form p n (s) = ±p„( 1 — s), 
and have zeros only on the critical line Re s = 1/2. 


1 Introduction 

Let Jb = 1, T\ = 1, A) = 2 and C\ = 1 and define the nth Fibonacci number, J~ n . 
and the nth Lucas number, C n , in the usual fashion [213, so that J ++2 = F n +1 + -F n , 
and the same recurrence for the Lucas numbers. We begin by examining the rela¬ 
tionship that exists between the Fibonacci and Lucas numbers, binomial coefficients 
and the Chebyshev polynomials of the first and second kinds, which we now define. 
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Definition (of Chebyshev polynomials). For 0 < 6 < ir, and n a non-negative 
integer, the Chebyshev polynomials (see for example mm) of the first and second 
kinds, T n (x) and U n {x), are defined by 

T n (cos9) = cos nO, U n (cos9 ) = ——-—(1.1) 

sm 6 

The related functions C n fix) and S n (x), are then defined to be C n fix) = 2T n fix/2), 
S n (x) = U n {x 1 2). 

Some key properties of the Chebyshev polynomials are given in Lemma 1.1, with 
equivalent Fibonacci and Lucas number expressions given in Lemma 1.2. 

LEMMA 1.1 (Chebyshev identity lemma). With T n (x ) the Chebyshev polynomial 
of the first kind, U n (x) the Chebyshev polynomial of the second kind, as defined in 
<o>, we have 



U n fix ) = 2 Tjfix), n odd, U n fix) = 2 Tjfix) — 1, n even, (1.4) 

j odd j even 



and 

T n+ ifix) = 2 xT n fix) - T n -i(x), 2T m {x)T n (x) = T m+n fix) + T\ m _ n \fix). (1.6) 

Proof. For proofs of the above identities we refer the reader to chapters one and two 

of pH. □ 

LEMMA 1.2 (Fibonacci identity lemma). With i 2 = — l, the Fibonacci and Lucas 
numbers can be expressed in terms of Chebyshev polynomials by 

77 _ Sn(i) r _ Cnfi) ^ _ n 1 o fi 

n+1 — . i L-n — -,,5 n — U, 1, 2, . . . , (1. 1) 

l n l n 

and in terms of binomial coefficients such that 



The binomial sum for the Fibonacci numbers is often referred to as the shallow 
diagonal sum of the Pascal triangle. The second sum can be viewed as a scaling of 
the terms in the previous sum, highlighting the interconnectedness that exists between 
the Fibonacci and Lucas sequences. 
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Proof. For proofs of these Chebyshev and Fibonacci identities, we refer the reader 
to pages 60-64 of m, where it is shown that one can deduce the binomial sums in 
(|1.8I) from the Chebyshev expressions in 4L7D - □ 

Definition (of cosine functions). For k € Z, let 4> m .k, ipmki fPnk and u rn k be defined 
such that 

<j) mk = 2cos (^pr), Vw = 2cos(^^r), 
pmk — 0m A; 2, — 0m k 2, 

and 

e(x ) = e 2mx = cos ( 2nx ) + i sin ( 2irx ), 

so that e{k/n ) is an nth root of unity. Then working (mod 2m + 1) the second 
Subscript, we have 0mj0m/c 0m(j-|-/c) ~b 0m(j—fc); as well as ^mffc+m+1) 0m/c; 
and in terms of the mth roots of unity 

4>mk = e ( 2m+l ) + e ( 2m+l) ’ V'mfc = e ( 4 ,^+ 2 ) + e ( 4m+2 ) ' 

Hence -0 2 2 = (1 + V5)/2 = 0, is the Golden Ratio, from which it can be seen that 

-021 = (l-\/5)/2 = 0 = 1/02 2- 

Definition (of vector convergents). We say that the sequence of vectors {v r }/i/ 1 in 
Q m converges to the limit v £ M m , if for any e > 0, there exists an r e € N, such that 
|v r — v| < e for r > r e , where | | denotes the standard Euclidean distance in R m . 

In this paper we identify the two (countably) infinite rational interlacing Fi¬ 
bonacci and Lucas sequences and J^/ 2,2 \ for r = 1,2,3,..., and their multi¬ 
dimensional analogues , with 1 < j < m. From ratios of these sequence terms 

we construct recursively (countably) infinite sequences of vectors \lR m € Q m , which 
converge (as per the above definition) to the irrational algebraic points = 

(0 m i,... ,0mm) £ IR m • These sequences of convergent vectors are a natural gener¬ 
alisation of the convergence of ratios of consecutive Fibonacci numbers Jv+i/J>, to 
0= (1 + >/5)/2, so that 

lim ~^~ r+1 = 0. 

1 >00 J- r 

To give a structural overview of this work, in Section 2 we use a rational recurrence 
relation to construct the two interlacing Fibonacci and Lucas sequences 


= £ 0 , -^1, £2, -F 3 , £4, -T. 5 , • • • = {£ 2r }^ 0 U {-^2r+l}“ 0 , 

r=l,2,3,... 

^ (2 ’ 2) =£i,-J- 2) £3,-J-4,£5,-J- 6> ... = {£2r+l}~ 0 U 1 , (L9) 

r=l,2,3,... 

(o 14 (2 9 ) 

as the numerator sequences of the respective rational sequences ’ ; and . 

The divisibility and continued fraction properties of these sequences are examined, 
and closed form expressions analogous to Binet’s formula are derived. The denom¬ 
inators of these rational sequence terms are powers of 5, thus indicating a 5-adic 
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structure. To achieve alignment of the denominator factors between the two se¬ 
quences, both sequences (as ordered above) commence with the term r = 1. 

Central to our results are the four families of polynomial functions P m (x), Q m (x), 
V m (x), and Q m (x), which are introduced in Section 3, with connecting trigono¬ 
metric, binomial, Fibonacci, Lucas and Chebyshev identities for these polynomials 
given in Theorem 1. The m-dimensional interlacing Fibonacci sequences of the title 
{,'Pr m ’ 3> }^L 1 , with 1 < j < m, along with the sister sequences are sub¬ 

sequently defined, and in Theorem 2 we obtain the generating functions for these 
sequences. As detailed in the Corollary to Theorem 2 and Lemma 3.2, this yields a 
number of ways to express (and such as 

m m / ■ )t \\ 

t= 1 (=1 C Sm Um+l^ 

With n = 2m + 1, the right-hand display allows for the geometric representation of 
the sequence terms as ratios of the diagonal lengths of an odd-sided regular n-gon 
inscribed in the unit circle (Theorem 5). 

The convergence properties of the pairwise ratios of these sequences are then con¬ 
sidered (Theorem 3), enabling the construction of the vector sequences € Q m , 

which converge to the limit point £ M m . Bounds for the remainder term and 
connections with simple continued fractions are briefly discussed in the two corol¬ 
laries. In Theorem 4 we show that the sequences (and are “rainbow 

sequences”, consisting of n times a renumbering of the Fleck Numbers (alternat¬ 
ing sums of binomial coefficients modulo n ) for r at negative integer values. The 
non-reduced numerators of are described in Theorem 6, where in 

particular for n = 2m + 1 = p an odd prime number, it is shown that 

jpL^J jr("M)) +0 ° = , so that G Z, V r £ Z. 

I J r=— oo v. J r=— oo 


In Theorem 7, we obtain sums of squares identities via the sequence term relations 




3 =1 


m, 0 ) 


3 = 1 


m d)V = 2 


2r+l i 


( 1 . 10 ) 


from which we deduce some combinatorial identities (not given in [6]), including 


E E (-ir +i+ “( 

j =1 \a =—oo 


2r+l 
r+j+a(2m+l) 


) = E <- 1 ) a ( 


4r+l 
2r+l+a(2m.-|-l) 


). ( 1 . 11 ) 


In Section 4 we examine in greater detail our Fibonacci, Lucas and Chebyshev 
polynomial functions, deriving orthogonality conditions, differential and recurrence 
relations, Mellin transforms with zeros only on the critical line lZe{s) = 1/2, a 
Christoffel-Darboux identity, and minimal polynomial relations. To conclude, in Sec¬ 
tion 5 we take a brief look at the sequences of matrix minors. 

Aside from the binomial Fleck number representation given in Theorem 4, as far 
as the authors are aware, the results contained herein are seemingly new and unpub¬ 
lished. Before deriving our results, we first reacquaint ourselves with the convergence 
properties of the Fibonacci and Lucas sequences. 
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Fibonacci Convergents 

The method known as Euclid’s algorithm or the highest common factor rule or the 
continued fraction rule is central to classical number theory. This algorithm produces 
all “good approximations” to a given real number a which can be rephrased in 
terms of 2 x 2 matrices such that, given a fraction a/c in its lowest terms so that 
hcf (a, c) = 1, find a matrix of integers 


a b 
c d 


with determinant 


a b 
c d 


ad — be = 1. 


As a consequence it means that much of number theory is fundamentally concerned 
with how pairs of integers behave, and so it is natural to study them accordingly. 

Taking a q = [a], the simple continued fraction algorithm produces a series of posi¬ 
tive integers ai, , ■ ■ ■, from which one obtains the convergents po/qo > Pi /qi , P 2 /Q 2 , • • •, 
such that po/qo = ao/1, and thereafter 


Pi 

qi 


1 


a o H-) 

a 1 


P2 

<?2 


1 

a o H-— 

ai -I- 

02 


Pn 

Qn 


— a 0 + 


ai + 


a 2 + 


1 

&n 


This algorithm allows the matrix representation 


far +1 Pr\ = fpr Pr-l\ (&r +1 A 
\^r+l Qr) \Qr Qr—l) \ 1 0/ 


{dr+lPr “1“ Pr—1 Pr 
T Qr—l Qr 


( 1 . 12 ) 


and by Dirichlet’s theorem for continued fractions (see m pi3i), the convergents 
have a remainder term satisfying 


Pr 

- a 

Qr 


1 1 

- < —• 

QrQr+1 q/ 


(1.13) 


As every real irrational number a has an infinite simple continued fraction expansion, 
one can think of R as being the completion of Q with respect to its Cauchy sequence 
limit points. 

The continued fraction expansion for the Golden Ratio (j) = (1 + \/5)/2 is [1; i], 
and so in terms of construction, it is the slowest converging simple continued fraction 
possible. From the expansion [1; 1] we deduce that the sequence of convergents obey 
the relation = 1 + with the first few terms in this sequence given below. 

1 2 3 5 8 13 21 34 55 
rr 2 ’3’ 5’ IF’ 13’ 21’ 34’"' 

Hence 


Pr r+2 r r+2 , 

— =- with hm -= 0, 

q r F r+ \ r^oo J r r+1 


Pr +2 Pr+1 T Pr 

Qr+2 qr +1 T qr 


(1.14) 


by (11.121) . concurring with the three term recurrence relation J - r +2 = ^>+1 + -T r . 
The corresponding recurrence equation is x 2 — x — 1 = 0, which has largest root 
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cj). As the sequence of Lucas numbers m obeys the same recurrence relation, with 
the initial values Cq = 2 and C\ = 1, we again find that the sequence of ratios of 
successive Lucas numbers converges to 0, and their reciprocals to 1/0 = 02i- 

In terms of 02 r , Binet’s closed form expression for the Fibonacci sequence, and 
the Lucas sequence analogue (see Theorems 5.6 and 5.8 p79 of [8]) equate to 


(_i \n 2 (_i \k+n 

Tn = (^2 - Ml) = E Mr, (1-15) 

and 

2 

C n = (-1 ) n (0? 2 + 0^) = E(-^r) n (1-16) 

r —1 

Subtracting the limit point of 0 = —0 22 from the sequence ratios J- n +\ /J~ n ] rewriting 
using the closed form (11.151) . and applying Dirichlet’s theorem (11.131) . we have that 
the Fibonacci convergent remainder terms satisfy 


' Fn+l 6 


v/5 



0 2n - (-l) n 


1 



Regarding the Lucas numbers for n > 3, and using (11.161) . we obtain the slightly 
weaker bound 


£n+i , 


V5 

1 

c * 

*-'71 


cj) 2n + (-l) n 

" r? 

L n-2 


Other variations on Binet’s formula, connecting the fifth roots of unity and the 
Fibonacci sequence, are 


(1 + (t> 2 k) n — n +1 + k €{ 1 , 2 }, (1-17) 

discussed by Grzymkowski and Witula in 0, and (-</>2 2 ) n = —02 2 -F n T -Fn— 1 : given 
by Vajda in [2Tj . 

From the relations given in (II.7|) . (11.161) and (11.171) one can deduce a multitude 
of identities, of which some of the better known are 

Fm-\-n -\-1 = Fm+lF rc+1 T F m F n , 

^n+l^n-l-^ = (-1)", 

Fm-\-n = n+1 F m —\F n) 

The two identities in (II.lip (proved in Theorem 7) are a variation on the right- 
hand display in () 1.18 1) . leading to representations of the sequence terms as 

a sum of the squares of m non-zero integers, and the sequence terms 2^1 T+i as a 
sum of the squares of 2m + 1 positive integers. The relations () 1.19 p and (|1.20l) high¬ 
light the interconnectedness that exists between these two sequences and also to the 
number 5. As a motivation for our general theories, we now describe an alterna¬ 
tive recurrence approach that yields two rational sequences, whose numerators are 
interlacing Fibonacci and Lucas numbers, and denominators powers of 5. 


n 


YA =F n F n+ l, 

(1.18) 

k=1 


Cn+lCn-1 - Cl = (—1)”“ 1 5, 

(1.19) 

F"m-\-n — §F m Fn-\-l Fm—1 F n . 

(1.20) 



































7 


2 Interlacing Fibonacci and Lucas Sequences 


Let = 2, 2 2 ’ 1 ' > = ~ 1) 2 ’ j) = 1) ^ 2 2 ’ 2) = — 1) and thereafter 

^ r + 2 ) = 3 € {1,2}, r = 0,1, 2 ,3,... 

so that in matrix notation 


/ «-(2,l) ®-(2,l) 

I ' 3 r+l 

*■(2,2) af (2)2) 

\ ^r+ 1 ^ 



1 

0 



(o . (o 9) 

where ; is the rth term of the first sequence and ; the rth term of the 

second sequence. Then both sequences satisfy a three-term recurrence relation, with 

recurrence equation x 2 + x + 1/5 = 0, whose roots are 


1-VE -1 , —1 — \/5 —(j) 

x = -=— =-= , and x = -=— = —=. 

2y/5 fy/E 2y/5 y/E 

(2 1) (2 2) 

LEMMA 2.1. The generating functions for the sequences TFr ’ ; and are 

given by 


2,1) . 5 (- x + 2 ) 
x 2, T 5x T 5 


3x 2 2x 3 7x 4 5x 5 18x 6 13x 7 47 x 8 34x 9 

2—x-|-1-1-1-h 

5 5 25 25 125 125 625 625 


Jf(2,2) . 


5x 


4x 2 3x 3 llx 4 8x 5 29x 6 21x 7 76x 8 55x 9 

. = \—x-\ -1-1-1-u 

x 2 + 5x + 5 5 5 25 25 125 125 625 625 


so that ^ r r ' 2 ' 1 ' > can be expressed in terms of J^ 2,2 ^ 


( 2 . 1 ) 


as 


«-(2,i) 


2^2) +^(2.2) 


r—1 • 


Proof. The result follows by applying the method of summation over the recurrence 
relation terms and rearranging. The relationship between the two sequences can then 
easily be seen by comparing the generating function structures. □ 


LEMMA 2.2. For r = 1, 2,3,... we have 



(2 1 ) (2 2 ) 

so that there exist closed form expressions for the sequences &)■ ’ ; and ’ ; in 
terms of tf 21 and fa 2 • 

COROLLARY. The sequence of ratios of consecutive terms with 

j £ {1,2}, approximates —cf/y/E with an accuracy 




r +1 


J?(2,l) 


+ 


V5 


1 


k2r-2 


+ 1 




r+1 , 

VE 


£■( 2 , 2 ) 

cJS 7* 



1 


4> 2r - 1 
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(2 2 ). 

Proof. We give the proof for ^ r v ' . Applying Cauchy’s Residue Theorem to the 

/o O'i 

generating function (12.11) of Pr ’ 7 yields 

f f(z) / 5z 

J F+ 1 d ~ = J z r+1 (z 2 + 5z + 5) d "’ 

where the contour contains the poles at z = 0 and z = , or equivalently 

z = —\/5^2i = 1^21 and z = \fh<f 22 = /a 22 . The residue at 0 gives the term TP^' 2 \ 
and as the sum of the residues is 0, one obtains 

jR2,2) = _£_ 5 

^ 21 (^ 21 -^ 22 ) ^ 22 (^ 22 -^ 21 ) 

Using 

M 21 - M 22 = \/5 and 0 2 2 = -—, 

021 

we deduce the desired closed form for P’r‘ 2,2 \ and similarly for Jpl 2 ' l \ 

The Corollary then follows from rearrangement of the closed form expressions of 
Lemma 2.2 in conjunction with —<t>/V 5 being the root of the recurrence equation 
with largest absolute value. □ 

LEMMA 2.3 (Numerator sequence lemma). Let jY^'^ be the non-reduced numer¬ 
ator of the rth term of the jth sequence , so that 

^( 2 ,i) = 5 L^J jK 2 j) 5 j € {i ; 2 }. 


with [J the floor function. Then for r = 1,2,3,..., we have 

P^r-l = T2r-2, ^ 2,1) = ^2r-l, and = C 2r - 1, = ^2r- 

(o I) (2 2 ).. 

Proof. Comparing the closed forms for TP? ’ ; and <Pf ’ ; given in Lemma 2.2, with 
Binet’s formula for the Fibonacci and Lucas numbers given in (|1.16l) . it then fol- 
lows that the non-reduced numerators of the two sequences ,Pf ’ ' and ,Pr ' ; are 
comprised of alternating Lucas and negative Fibonacci numbers (by non-reduced 
we mean that any common factors between numerator and denominator, such as in 
5/25 in the sequence term P r ^ ’ , are not cancelled). □ 


As described in the Corollary to Lemma 2.2, the ratios of consecutive sequence 
terms yields two sequences of convergents with common limit point —<f>/y/§. The 
initial terms of these two convergent sequences are given below. 


^(24) 

^r+1 

1 

3 

2 

7 

5 

18 

13 

47 

34 

123 

89 

7 * 

2’ 

5’ 

3’ 

10’ 

7’ 

25’ 

18’ 

65’ 

47’ 

170’ 

123’ 

«-(2,2) 
^r+ 1 

= 1, 

4 

3 

11 

8 

29 

21 

76 

55 

199 

144 

jS-(2,2) 

5’ 

4’ 

15’ 

11’ 

o | 

29’ 

105’ 

76’ 

275’ 

199’ 


LEMMA 2.4. Allowing negative integers ai, recurring integer continued fraction 
expansions for the two sequences TP}■ 2 ' 1 ' 1 and TP^'^ are given by 

^ = [0; 1,2,3, —i, 5], ^ = [1; —5, i]. 
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(2 i') 

Proof. For JFf ’ we assume that [ 00 , 01 , 02 , 03 ] = [0,1,2,3], and then apply the 
continued algorithm to deduce that 04 = ae = —1, and 05 = 07 = 5, so that the 
expansion recurs from that point onwards. The continued fraction expansion for 
&r ’ 7 can be deduced similarly. □ 

Definition (of divisibility sequences). Let 01 , 02 , 03 ,... be a sequence of integers 
satisfying the divisibility property s \ t => a s \ at, with s and t positive integers. 
Then we say that the sequence of integers {a r }^ =1 is a divisibility sequence. Moreover, 
if we have the stronger condition 

hcf(s, r) = d =>• hcf(a s , a r ) = o^, 

then we say that it is a strong divisibility sequence. 

LEMMA 2.5. Let s and t be positive integers with s \ t, so that t/s is an integer. 
Then M^ 2 ’ 2 ^ \ JY^ 2 ' 2 ^, and the numerator sequence {Mr ( ' 2 ’ 2 ' ) }£L 1 is a divisibility 
sequence. Moreover, we have that M^ 2,2 '* = JY^^ jYr 2 ' 2 \ and if s/d and t/d are 
both odd integers, then hcf(s,t ) = d hcf^PY^ 2 ^, 

As a preamble to proving this lemma we now introduce the Gauss hypergeometric 
functions, as well as the Fibonacci and Lucas polynomials, which underpin many of 
the results discussed later on. 


Definition. For integers r, k , with k > 0 let the rising factorial, falling factorial and 
hypergeometric functions be respectively defined in the usual manner such that 


r k = r(r + 1 )... (r + k — 1 ), r- = r(r — 1 )... (r — k + 1 ), 


and 


( 2-1 , , ( 2-77 

61, . . . , 677 


* = 




where 


tk = = 


a\ ... a k „z k 


k >0 b l ■ ■ ■ b n kl 

with none of the hi zero or a negative integer (to avoid division by zero). 


Definition (of Fibonacci and Lucas polynomials). Let the Fibonacci polynomials 
J- n (x), be defined by the recurrence relation 


F n+ i(x) = xF n (x) +F n -i(x), with Fi{x) = 1 , F 2 (x) = x, ( 2 . 2 ) 


or equivalently by the explicit sum formula 

= x"- 1 2 Fi f 


[(*»-1)/2] 
Fn{x) = ^ 

3=0 


n-j-1 

j 


n n 

—, 1 -; 1 — n; — 

2 2 



(2.3) 

Similarly, define the Lucas polynomials C n {x) by the recurrence relation in (12.21) . 
but with initial values C\{x) = x, C 2 (x) = x 2 + 2 , or equivalently by the explicit 
sum formula 


C n {x) 


\n/2\ 

E 


j =0 


n 

n-j 




= x n 2 F 1 




(2.4) 
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LEMMA 2.6 (Factorisation and divisibility lemma). The Fibonacci and Lucas 
polynomials and numbers have the following factorisation and divisibility properties. 

(1) With i 2 = —1, the roots of the Fibonacci and Lucas polynomials can be ex¬ 
pressed in terms of i multiplied by the cosine of rational multiples of it, so that 
the polynomials can be factorised as 



(2) The Fibonacci and Lucas polynomials have the divisibility properties 

F n (x) | F m (x) &n\m, F n (u p -i (\/5/2^ = T np /T p , (2.7) 

and 

C n (x) | C m (x) m = (2k + l)n, for some integer k. (2-8) 

(3) The Fibonacci sequence is obtained by setting x = 1, so that T n = F n (l), and 
similarly for the Lucas sequence with C n = C n ( 1). Taking x € {2,3,4,... }, 
then produces one possible definition of higher dimensional Fibonacci and Lu¬ 
cas sequences which obey the divisibility properties as stated. 

(4) Two important divisibility characteristics of the Fibonacci sequence are that 
for p a prime, p divides ^ s y with (j^j the Legendre symbol, and that 

hcf(s,r ) = d =>■ hcf(F s ,F r ) = T& so that it is a strong divisibility sequence. 

(5) With the extra constraint that s/d and t/d are both odd integers, an analogous 
divisibility sequence result holds for the Lucas numbers whereby hcf(s , r ) = 
d =>- hcf(C s ,Ct) = Cd- It follows that if s/d is an odd integer, then Cd divides 
C s . The Lucas numbers also have the factoring property F-z r = F r C r , so that 
C r | r- 

Proof. For proofs of these divisibility properties see T. Koshy [8], pl96-214, and 
p451-479. □ 

3 Higher-Dimensional Interlacing Fibonacci Sequences 

In order that we may generalise our previous results for FI 2 ’ 1 ' 1 and Fr 2 ' 2 \ to 
with 1 < j < m, we now introduce families of polynomials related to the Fibonacci 
and Lucas polynomials, that are central to our theories. 
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Definition (of generating function polynomials). For positive integer m, we define 
P m (x), Q m (x), V rn (x), Qm(x ) and V m (x), to be the polynomials of degree m given 
by 


/ \ \ ^ 2 m + 1 (m k\ u „ , , \ - 

P m(x) = Y ( 2k Jx , Q m (x) = 


k =0 


m f m + k — 1 


k =o 


Pm{x) = Y 


k =0 


m + k 
2k 


Qm(x ) = Y 


k =0 


k \ 2k — 1 


m + k + 1 
2/c + 1 


x 


x k , (3.1) 


(3.2) 


and 


7T2 / r L* 771 ”1 \ 

u.w = E(- 1 )” +[!+?l ( [5 t Tl ) i ‘- 

fc =0 ' ' 


where the identity 


lim 


j (j + k- 1 


= 2 , 


(3.3) 


(3.4) 


fc->o k V 2/c - 1 
ensures that Q m (x) is well defined. 

We label the roots of P m (x), ordered in terms of increasing absolute value, by 
■ ■ • ,17mm, and similarly i/ ml , v m2 ,. ■ ■, v mm the ordered nn roots of Q m (x), 
so that we may write 


Pm{x) — J^(x Qm(x') — |^(x J'mi), 


(3.5) 


2—1 


2—1 


where for i < j, we have \y m i\ < |/J mj |, and | Vm i | 5; | k'm j | • 


In Theorem 1 we show that the above definitions for y m i and v m agree with those 
previously given in the definition following Lemma 1.2. This leads to simple identities 
for the polynomials P m (x). Q m (x), V m (x), Q m (x) and V m (x), in terms of Chebyshev 
S m (x ) and C m (x ) polynomials, as well as Fibonacci and Lucas polynomials. 


THEOREM 1. The roots of the equations P m (x) = 0 and Q m (x) = 0 are real, 
simple, negative, contained within the interval [—4,0], and with the above definitions 
for p m k and v m k, 1 < k < m, we have 


Tm k — fim k 2 — 2 COS 


\2m+l) 


so that P m (x) = 


Vrn k — firn k 2 — 2 COS 


( ~ X ) \ _ 9 
V 2 m+ 1 ) 

(3.6) 


m / / 2nk \ \ 1 1 

JJ ( x + 2 — 2 cos ( 2 m + 1 ) J = T7r^ 2m+1 = 7/x 2m + 2 (v^) + p 2 m{Vx)) , 

k =1 V V / / V V 


(3.7) 


= U 2m 1 + ^ = S- 2 m (Vx + 4) = S 2m (2 cos y ), with x = 2 cos 2 y - 2, 


(3.8) 


Qm(x) = ( x + 2 — 2 COS 


k =1 


ir(2k — 1) 
2m 


— P-2in (y/x'j — P 2m -\.\( y \/x)- i rP 2 m — l{.''J~x) 

(3.9) 
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= 2 T 2m (^J 1 + ^ = C 2m (\/x + 4) = C 2m (2 cos y), with x = 2 cos 2y - 2, 

(3.10) 


and 


V m (x) = JJ ^x + 2 + 2cos ^ 2m + l )) = ' F ‘ 2m+l (^) ’ ( 3 - n ) 


= G 2m ( —^ ] = Sim W~x) = Sim (2z cos y), with x = 2 cos 2y + 2, (3.12) 



Qm{x) — /—•F 2 m+2 (\/®) > 

V® 


and 


y m (x) = P m (x - 2) = ( x - 2 cos 


k =1 


2irk 
2m + 1 


(3.13) 

(3.14) 


— G 2m 


l + x/2 


= (V* + 2) = S 2m (2 cos y), with x = 2 cos 2y. 

(3.15) 


COROLLARY. the also have the relations 


xPm— 

i(*) 

— Qm{x) Qm-l(^)? 

Qm{x) — 

P m (x) 

— 

Pm— \ (x ), 

(3.16) 

xVm- 

l(*) 

= Qm{x) Qm—l{x)-> 

Qm{x) = 

V m (x) 

- 

Vm-l(x). 

(3.17) 



Vm(x) = (-ly 

m Pm{-X~ 

4), 



(3.18) 

xP m (-x 2 ) 

= (- 

-l) m—1 Q 2 m+l( x 2), 

Qm ( X 

! ) = (- 

d ; 

m Q2m{-X 

-2), (3.19) 

and 









Qm( 

J x P mi [x)V mi (x) 

X) = < 

if m 

= 2m i 

+ 

1 is odd, 

(3.20) 



^ Qm\-\-l (x) Qmi (^) 

if m 

= 2m i 

is 

even, 


along with the integral identity 


(2m+ 1) f V m (t 2 )dt = y/xP m (x) => (2m + l)V m (x) = P m (x)+ 2xP’ m (x). (3.21) 
do 

Proof of Theorem 1. For A; = 1,2,..., m, the function cos ^ 2 m+i ) a decreasing 
function of k. It follows that cj) m k is a decreasing function of k, and so in terms of 
absolute values we have, /j m \ < y m2 < ... < y mm . A similar argument holds for 
the 1 < k < m, and hence (13.51) . 

Although it follows from (13.61) that the roots of P m (x) are simple and lie in the 
interval [—4,0], we demonstrate this by two other methods in order to highlight the 
links that exist between P m (x ) and the Legendre and Chebyshev functions . 

Method 1. By using manipulations with Pochhannner symbols which we omit, 
P m (x) may be written in terms of the Gauss hypergeometric function 2 Fi. Letting 
P™(x) denote the associated Legendre function, in turn P m (x) may be expressed as 


(2m + l)y / 7T (x + 4) 1 / 1 „_ 1/2 / 
l( 1 2 (-x)V4 l 


1 + f 


(3.22) 
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_^ /2 

The functions P m (z) are orthogonal on the interval [—1,1]. With x = 2(z — 1), 
it follows from a standard result in the theory of orthogonal polynomials [20] that 
the zeros of Pm l ^ 2 (l + f) are contained in [—4,0] and simple and hence for P m (x) 


too. 


□ 


Method 2 It is known that (e.g., [12| . p. 64) 

2 sin[(t' + l/2)y)] 


P~ x ' 2 (cos 99) = 


7T sin ip (u + 1/2) 


□ 


Then with p replaced by cos -1 ip and x = 2(ip — 1) it again follows that the zeros of 
P m (x) are in [—4,0] and simple. 

Remark. It is then possible to write 

where U u _ 1 / 2 is the Chebyshev function of the second kind. 

To see (|3.14l) and the first part of (13.71) . we have 


k =1 


( x — 2 cos 


2nk 
2n + 1 


2 n 


k =1 


X 


= Un hr + Un-l hr = U 2n 


— COS 


1 + x/2 


2tt k 
2n + 1 


1 + 2 Tk 


k =1 


= Pn(x- 2) = v n (x), (3.23) 


where we have used the relations (|1.3I) . (11.41) and (jl.5D . of Lemma 2.1. 
Hence 


k =1 


j J f x + 2 — 2 cos 


27 rk 
2n + 1 


l + 2^2 T k (- + 1 


k =1 


— Un ( 77 + 1 ) + U n —1 ( — + 1 ) — Pin A/ 1 + 7 — Pn{x), 


(3.24) 


and the result follows. Similar arguments can be used to obtain the first statements 
of (13.91) and (13.111) . 

The latter statements in (13.71) . (13.91) . (13.111) . and that of (|3.13p . concerning ex¬ 
pressions for P m (x), Q m (x),V m (x) and Q m (x) in terms of Fibonacci and Lucas poly¬ 
nomials, can be deduced via binomial relations as follows. We have 


1 1 m 

Vi 

m 

-£ 


2m + 1 f 2m + 1 — j 


2m + 1 — j 


J 


(VS) 


2771+1 — 2 ^' 


3=0 


2m + 1 — j 


2m + 1 (2m + 1 — j , 

x J = 


£ 

3 =0 


2m + 1 fm + j 


2 j + 1 \m — j 


2m + 1 fm + j 


E Z.IU -|- 1 

2j + l 


j=o 




2j 

2m f 2m — j 


^ = P m (x), 


3 =0 


2m — j 


(VS) 


2m—2 j 


and 
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y 2m (2 nn-j\m-i 
y 2 m - j V j ) 


E m 

j 

j =0 J 


m f m + j — 1 

2 j-l 


Qm(®) ■ 


The Fibonacci polynomial identities follow similarly. The respective generating func¬ 
tions for the Lucas and Fibonacci polynomials (see page 447 of 0 ) are 


G L (x,t ) = Y2 £n(x)t n 

71—0 


1 + t 2 

1 — t 2 — tx ’ 


<MM) = ^-^n(x)f n = 2 

71—0 


so that tGi,{x,t) = (1 + Hence we have C n {x) = J r n +i(^) + .F n _i(x), 

giving 

Pm{x) — 77 [•^2m-l-2(V^c) T •^2m(\/^)] j Qm{x) — •^2m+l(v^) T -Fltn— 1 (’s/x): 

V ® 


and so the ordinary generating functions of P m (x) and Q m (x) can be written as 

oo oo 

Gp(x,t ) = ^P n (x)t" = i_(2+*)t+P ’ Gq( x G) = y; Qn{x)t n = i_(2+t)t+p ■ 


71—0 


71—0 


(3.25) 

Combining the above relations with the identity P 2m {x) = P rn (x)C m (x), we obtain 

Pm{x)Pm{x) — ^^2771+1 (V^)^2m-|-1 (\/®) — ^-^4m+2(\/^) — Q2m{x ), 

VI v x 

Qm-|-l(l)2m(l) — ~^P-2m+2 (’^/x)J~2m+2 (y/x) — ^J r 4 m -|-4(\/x) — Q2m+1 (®) • 

Vi VI 


The Chebyshev identities in (13.81) . (13.101) . (13.121) and (13.151) can be derived directly 
from the definition (tm with further connections to the Chebyshev polynomials 
established via the expression for P m {x) given by 


Pm(x) = T m (l + x/2) + yfl + 4/xsinh[2m csch l (2/yfx)\. 


The expressions (13.161) and (I3.17P of the Corollary, follow from direct manipulation of 
the binomial coefficient polynomials (see Lemma 2.2 of [4]) given in the definitions for 
P m (x ), Q m (x), Vm(x), and Q m (x). To see (I3.19p . substituting — (x+ 2 ) in the product 
formula for P m (x) in (|3.7I) and comparing with (13. 14p . gives us P m (x — 2 ) = V m (x), 
and writing 


Q2m+l( (x p 2 )) 


2771+1 

n 


k=l 


x — 2 cos 


7 r( 2 fc - 1 ) \\ 

4m. + 2 ) J 



= (~ x )Y{(x 2 -2-4> mk ) = (-l) m 1 x Yl (~x 2 + 2-<f> mk ) = (-l) m 1 xP m (-x 2 ), 


k =1 k =1 
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we obtain the first expression in (13.191) . Similar arguments produce (13.181) . the second 
expression in (13.191) . and (J3.20I) . Considering (13.211) we have 



m 


*E 


r+‘) ( “ 


2k + 1 


y/x 


t =0 


Pm ("0 ■ 


and differentiating we obtain the final display. 


□ 


Definition (of m-dimensional interlacing Fibonacci sequences). Let the matrices of 
binomial coefficients B 0 ^ and B even be defined such that 

-Bodd = ( bi,j)mxm, with b h3 = (-1)* +J_1 , 

and 

Beven — (^ij)mxmj with bjj — ( 1) ^ 

and let the recurrence matrix R m be given by 



/ -hi 1 0 ... 0 \ 

-h 2 0 1 ... 0 


—h m -1 0 0 0 1 

V -hm 0 0 0 0 / 


with 


1 / rn + k\ 

k = 2fc + 1 \ 2k )' 


(3.26) 


Also let the sequences of m x m matrices M 0 (m,k) and M e (m,k ) be respectively 
defined such that 


M 0 (m, r) = (2m + 1 )B odd R™ +r , and M e (m, r ) = (2m + 1 )B even R^ +r , 

so that each row of M 0 (m, 0) or M e (m, 0) corresponds to a list of consecutive se¬ 
quence values. The matrices M a (m, 0) or M e (m , 0) are then taken as the initial value 
matrices for the two sets of m interlacing Fibonacci sequences that they generate, 
under repeated multiplication on the right, by the recurrence matrix R m . 

We denote by , 1 < j < m, the rth term in the jth rational interlacing 

Fibonacci sequence of dimension m, and by 1 < j < m, the type jth 

rational interlacing Fibonacci sequence of dimension m. The case m = 2, then cor- 
responds to the interlacing Fibonacci and Lucas sequences ; and J^ r v ’ , defined 
in Section 2, and for the case m = 1, we simply have 

^ (11) = (if”) ’ ^ = 1,2,3,... . 

It follows that row j of M 0 (m,r ) contains the jth m-dimensional rational interlac¬ 
ing Fibonacci sequence terms • ■ • > and row j of M e (m,r) 

contains the type jth m-dimensional rational interlacing Fibonacci sequence 
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terms ( ^r+m > ^r+rn- 1; • ■ ■ ■ i where both sequences satisfy the recurrence re¬ 

lation (given here in terms of 


'^r+m 


1/m + l\ ( mj ) ( m + 2\ y(mj-) _ 1 I'm + m\ ar ( m ,r) 

3 \ 2 J^r+m-1 5 ^ 4 J^r+m- 2 ■■■ 2 m + 1 V 2m y r+m ■ 


For example, when m = 

5 and 

r = 

2, we 

have 









( "I 

3 

-10 

35 

-126 \ 


/ 

-5 

1 

0 

0 

0 \ 


0 

-1 

5 

-21 

84 



-7 

0 

1 

0 

0 

M 0 (5,2) = (2m+ 1) 

0 

0 

-1 

7 

-36 



-4 

0 

0 

1 

0 


0 

0 

0 

-1 

9 



-1 

0 

0 

0 

1 


V 0 

0 

0 

0 

-1 / 


V 

I 

II 

0 

0 

0 

0 / 


7 


( jf - 1 ) 

^(54) 

^6 

... jf ’^ \ 

/ 

®-(5,2) 

^7 

*■(5,2) 

^(5,2) 

... ^3 


®-(5,3) 

^7 

df(5,3) 

^(5,3) 

... ^3 

= 

^•(5,4) 

^7 


... ^-3 


^(5,5) 

jf(5,5) 

^■(5,5) 

... ^3 / 

V 


10744 

11 

28817 

37V34 

341)69 

20&2 

11 


3415 

11 

99 

-32 

11 ^ 

9156 

H 

265 

-85 

28 

11982 

llV02 

6^5 

11 

346 

-110 

35 

317 

-100 

31 

188 

-59 

18 ) 


LEMMA 3.1. The product of the eigenvalues of the recurrence matrix R m (and 
so its determinant) is given by (—l) m /(2 m + 1), and the sum of the eigenvalues 
by —hi = — m(m + l)/6. The binomial matrices satisfy Det(B oc jd) = (— l) m and 
Det(B even ) = l, so that the two sequences of determinants of the matrices M 0 (m,k ) 
and M e (m , k), for k = 0,1, 2,3,.. consists of terms of the form ± • 

The inverse matrix of R m is given by 



( 0 

0 

0 

0 

— (2m + 1) 

\ 


( 0 

0 

0 

0 

~9m 'l 


1 

0 

0 

0 

— (2m + l)hi 



1 

0 

0 

0 

9m— 1 

TD— 1_ 

lx "m 

0 

1 

0 

0 

— (2m + 1 )h 2 


= 

0 

1 

0 

0 

9m—2 


\ 0 0 0 1 —(2m + l)h m -i / V 0 0 0 1 ~9i 


where 

2m + 1 / 2m — k 

9k ~ 2m + l-2k V k 


so that gk = (2 m + 1 )h m -fc- 


(3.27) 


Denoting by R m (x) and R rr ) (x) the respective characteristic polynomials of R m and 
Rif), we have 


Rm(x) — ^ " h m -j x^, H m (*^) — (2m -(-1) ^ ( hj xj — E m (a;), 

3 =0 3 =0 


where the eigenvalues of the inverse recurrence matrix Rf) are the roots of the 
polynomial P m (x) and the eigenvalues of the recurrence matrix R m are the roots 
of the polynomial P“ 1 (x) (say). By \3. 6)1 of Theorem 1, the roots of P~ l (x) are 
given in descending order by p m \, Mm2• • ■ ■ ■ Mmm• Hence, with a suitable choice of 
(algebraic number) coefficients we can write 


&(rn,j) _ (m,j) -r 1 J m J) ..-r , 0 (w) -r 

— u 1 Mm jt «2 Mm 2 ' ■ ■ ■ “m. Mm rr 


(3.28) 
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Moreover, if 2m + 1 is an odd prime number, then R™ contains integers above the 
principal diagonal and rational numbers with denominator (2m+l) in its lower trian¬ 
gular matrix part. Consequently, the initial value matrices M 0 (m, 0) and M e (m, 0) 
for the respective sequences and are integer matrices if and only if 

2m + 1 is an odd prime number. 


Proof. For a general m x m recurrence matrix K m of the form 


K — 


— ai 

1 

0 

0 ... 

0 

-02 

0 

1 

0 

0 

Am-1 

0 

0 

0 

1 


0 

0 

0 ... 

0 / 


(3.29) 


evaluating the determinant along the first column yields \K m \ = (—l) m a m , and 
similarly we obtain the characteristic polynomial K m (x) = 


\K m — xI m \ = (-l) m (a m +a m - 1 x + .. . + aix m 1 +x m ) = (-l) m ^ a m -j x j . (3.30) 

j =o 

Hence if a m 0, then the inverse recurrence matrix R~ 3 exists, and is given by 

/ 0 0 ... 0 --L \ 

Uj m 

1 0 ... 0 —2L 

K - 1 = 0 1 ... 0 —Si 


o o ... 1 - 


&m — 1 


with characteristic polynomial 
(- 1 ) 


Km{x) = -—(l + a\x + ... + a m -ix m 1 + a m x m ) = -—V ajX j , (3.31) 
a m a m ^ J 

j=0 

and the expressions for R m , i?” 1 , R m (x ), and Rf n l (x) follow. 

If the m eigenvalues of the recurrence matrix K m , denoted by A m i, A m 2 , ■■■, A m m , 
are non-zero, real, algebraic, distinct, and listed in descending order in terms of ab¬ 
solute value, then the sequences 


YnP — (l) \nji ^mji ^mji ■ ■ ■) > j ~ ^- 1^1 • m, 


form a basis for the solution space of all possible sequences satisfying the recurrence, 
for any possible initial values. Similarly the ordered eigenvalues A~’ m ,..., A“\ of the 
recurrence matrix Riff, form a basis for the solution space of all possible sequences 
satisfying the inverse recurrence relation. It follows that each sequence term gener¬ 
ated by the recurrence relation can be expressed as a linear combination of powers 
of the eigenvalues of the recurrence matrix. With regard to our sequences 
and recurrence matrix R rn , we have A m j = 1 /fa m i, and hence (13.281) . 
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The statements concerning the entries of R™ and that M 0 (m, 0) or M e (m, 0) are 
integer matrices when 2m + 1 is a prime number follow from the property 

2 m + 1 (m + k\ 

2k + 1 \ 2k ) € ' 

□ 


THEOREM 2. For m a positive integer, and 0 < j < m — 1, we have 

o / • , , , i \ j 




k=0 


j + k + 1 \ 

2fc + i y r - fc ’ 


m,m—j) _ _ Sf' ( .J + k \ 

2/c / ^ 


E 

k=0 


(3.32) 

so that each of the terms in the sequence can be expressed as a binomial 

coefficient linear combination of (m + 1 — j) terms from the sequence 
where we note that = V r € Z. 

T/ie generating functions for ^ and ^ are given by 




m ~j)x r = 


r=0 


(2m + l)Qj(x) 
Pm(x) 


Y.F 

r =0 


m ,m-j) x r — — 


(2m + l)'Pj(x) 
P,n(x) 


(3.33) 

so that the sum of the numerator coefficients of each generating function is a Fi¬ 
bonacci number. 

Inverting the expressions for and&r™’^ in ()3.32|) . we obtain 


J 1 2j -1 fj + k~l\ {m , 1) 

2fe ’ - J - m ’ 


fN) _ V —_3: 

r ^2k + l 


k =0 


and 


r 


fc =0 


j {j p k l\ ~(m, 1) , , • , . 

jfe( 2*-l F-* ’ iSySm-i. 


(3.34) 


(3.35) 


COROLLARY, /n terms of the roots Umii ■ ■ ■ ,hmm of the polynomial equation 
P m (x) = 0, we have 


9m 4 - 1 

&r(m,m) _ _c#(m,m) _ _ ~ l ~ 1 _ 


-(2m+l)X 2 


k =1 


I COS 


2nk 
2m + 1 


2 II 2 




i cos 


2nk 
2m + 1 


= ,, 1 - r _L _i_ // 1 =X^ I 9 

^r rml W • • • W P'mm / . ( ^ 

fc=l 


COS 


2irk 
2m + 1 


— 2 cos 


1—r 


(3.36) 
2vr j 


— 2 


2m + 1 

(3.37) 

(3.38) 


-l 


j - 1 


^ (raj) = X (^t) 1 r <Ani) = X (Vmt) 1 r 0 0“™* ~ M(j-l)fc) 


t=l 


t=l 


fc=l 


J-l 


XI r II ~ = X (^i) 1 r S '2j-2 ( 2 


t=l 


fe=l 


i=l 


I COS 


TTt 


2m + 1 


(3.39) 
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= (-ir 1 E 


2 sin 


2m+l j 


t =i ( 2 sin 


/ 7T t 

y 2m+l 


2r—1 


E^ m t ) 1_r y*- 1(0 mt) ? 


t=l 


and 


^ = E (/^t) 1_r Qi (fet) = E ^rnt) 1 - r n ^ 


m t 


Vjk) 


t =1 


t=l 


fc=l 


x> 

t=l 


m t J C 2 j 


2 COS 


nt 


2m + 1 


t-ir'E 


2cos (sSt) 

2r^2' 


t=i (2 sin 


/ TVt 

y 2m+l 


(3.40) 


Proof of Theorem 2. The denominator polynomial of 2 m.+\ Pm(x) in the generating 
functions for (and so ^”'•"4) i s a direct consequence of 2m \i Pmfx) being the 

recurrence polynomial for . The simple numerator follows from the starting 

vector for the initial value matrix consisting of (0, 0,0,, —1). 

The two identities in (13.321) follow directly from (6.8) and (6.9) of Lemma 6.2 
in m- Applying (13.321) to the generating function polynomial for thus es¬ 
tablishes the numerator polynomials in the generating functions of ^ and 

in ([3733]) 

With the binomial matrices of initial conditions M 0 (m,r) and M e (m,r), so de¬ 
fined, it is possible to invert the identities in (13.32[) using a binomial convolution to 
obtain (I3.34|) and (13.351) . 

To see (13.361) . we know that P m (x) factors as P m {%) = n^LiO^ — Tmk)- Our 
contour integrals from the generating function have the form 

1 r (2m + l)z 

f z^P m (z) dZ - 

The contour encloses the origin and at least the interval (—4, 0] of the negative axis, 
in order to contain all of the simple poles of P m (z) and the higher order pole at the 
origin. 

Assuming the partial fractional decomposition 

1 


Pm{x) 


E 

k =1 


Ck 


X Hm k 


we now show that Ck = 1 /P' m {Hmk), where the condition that P m (x) has distinct 
roots implies that P' m {^ m k) y 0. Using the form with lowest common denominator 
P m (x), we have 


P m {x) 


E 

k =1 


Ck _Yhk=l c kWj=l,j^ki x Pmj) 

X hmk Pm(x) 


Then 1 = TZ=i<*n?=i J¥ *(* — Hmj ), and, evaluating at fi mn , 1 < n < m, gives 
f c n DIj —1 j^ n {hmn hmj) — c n P m (^fi mn ) . Hence c n — \/P rn (^fi mn ). 

We have determined that 


E 


Ck 


Pm(pc) ^ l^mk Jq—i ^ ^"mk) Pm^mk) 


E 
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wherein P^{n mk ) = UjUj^ki^mk ~ Vmj)- Then 

1 i 2m + l d _ 1 /2m + ly 1 1 ^ 

2m J z r+1 P m (z ) * 2ni J z r+1 ^ (z - /j, mk ) Pm(Vmk) 

and the residue at 2 = /i mk is given by 

(2 m + 1) 1 (2 m + 1) 1 

^mk P m{Hm k) Hj=l,j^k(^rnk ~ ^mj) 

The pole at the origin gives the term generally. Using 27ri times the sum of 

all residues gives 




2m + 1 


^mk n tfkiv m k f^m j ) 


= 0, 


r = 1,2,3,..., 


(3.41) 


and hence the result and (13.371) . 

The identity + P m 2 + • • • + I- 1 mm i n H13.38I) follows similarly to 

(13.361) . Combining these results with (13.341) and (|3.35l) then gives 


= f y ~ 1 (i + k 

r ^^2k + l\ 21 


- 1 


t =1 k =0 


2k 


2 cos 


( 2nt \ 
\ 2m + 1 ) 


-2 


k+l—r 


1 < j < m. 



m m j— 1 

= /* ] ik-mt) P j — 1 (pmt) = /* ] (l^mt) | (f^mt k-(j—l)k) > 

t= 1 t= 1 k =1 


and 


(g(™,j) _ L 

r ^ k 

k =0 


j + k - 1\ 
. 2fc -1 y 


^r—k ’ 


1 < j < m — 1, 


2 j / j + k — 1 
2k \ 2k-l 

t =1 fc=0 v 


2 cos 


/ 2vrt \ 

V 2m + 1 / 



fc+1—r 


m 


-E 


( 2 cos ( 


2-Kt 

V 2m + 1 


)- 2 


1—r J 


2 j_fj + k 


- 1 

^ 2k\ 2k - 1 

fc=0 v 



/ 27Tt \ 
1 2m + 1 / 



^2 cos 


2 cos 


2irt 

2m + 1 


2nt 

2m + 1 




k 


m m j 

^ ' (gmf) Qj ( k'mt) ^ • 

t=l i=l fc=l 

We have thus established (13.391) and (13.401) . 

□ 
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Remark (to Theorem 2). For 1 < j < m, the sequence terms and 

naturally occur in matrix powers of particular circulant matrices as described in 
m- In this setting there exists an additional sequence which appears in 

the leading diagonal of the powers of this circulant matrix, thus bringing the total 
number of sequences to n = 2m + 1. This sequence satisfies the relation 

m 

= 2 jK™d) = -2^2 (3.42) 

k =1 

so it also obeys the recurrence relation R m , and we can write 




+ k ~A 

^ k V 2k - l J 

k =0 v ' 


A m P) 

*r—k ’ 


1 < J < 777. — 1. 


(3.43) 


By Lemma 6.1 of m, (see also Lemma 3.3 of mi) these type of recurrence relations 
can be written as half-weighted determinants. For example has the form 






(T) 

(T) 

i 

0 

0 

0 

... 0 

hi 

1 

0 

0 

... 0 

(T) 

h 2 

hi 

1 

0 

... 0 

fm+r—1\ 

V 2r-l ) 

h r —i 

hr-2 

h r -3 

hf — 4. 

... 1 

(m+r\ 
\2r+l/ 

h r 

h r — i 

h r -2 

h r -3 

... hi 


Hence for r € N, the sequence terms 1 < j <,m, are rational numbers 

corresponding to a rational multiple of either the r or the (r+1 (-dimensional volume 
of the simplex with vertex coordinates given by the row entries of the determinant. 


LEMMA 3.2 (Ratio lemma). For j > 0, and 1 < t < m, we have 


hmt Pj(hmt) ~ hm(j+l)t hmjti hmt Qj{hmt) ~ hm (j+l)t 2/7 m jt + 

(3.44) 

and 


hm (m-2j -\-1 )t hm(m—2j)t _ 4 > m(m—2j-\-l)t ( t > m(m—2j)t 

hm (m—j+\)t hm (m—j)t Pm (rn—j-\- ] )t Pm (m—j)t 


fimjti 


hm(2j—m)t hm(2j—m—l)t Pm. (2j—rn)t 0m (2j—m— 1)7 
hm(rn—j+l)t hm (m—j)t Pm (rn-j-\-l)t Pm(rn—j)t 


Omjti 


Moreover we have 


1 < 3 < [f ], 

(3.45) 
[f] + l <j<m. 

(3.46) 


j - 1 

hmt | (hmt h(j—l)k) = hmjt ~ hm(j—l)t = fimjt ~ 0m(j — 1)7) (3-47) 

k =1 


j 

hmt | {hmt Vj k) = (f> m (j-|-l)i — 20 m jj + 0m (j— 1)7) 

k =1 
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so that 


^ (raj) = (Vmt) '' {</>m jt - 0m (j-l)t) 7 
t= 1 
m 

^r = 'y y (fbnt) (0m (j+l)t ~ 2 (f>mjt 4" ■ 


(3.48) 


(3.49) 


t=i 


Definition (of convergent vector sequences in Q m ). We define the convergent vector 
sequences . and their limit points such that 3>( m ) = (0 ml , 0 m2 , ■ ■ ■, <A m m) • 


and = 






-&r 


-&r 


— c^r* 


& 


( m,m ) ’ ‘ 


g-(m,m-r+ 1) ’ " ' ^-(m,[m/2]+l) ’ ' 

r r 


1) 




(m,l) 


where from the definition of , we can replace by in the above 

expression for . without affecting the ratios. 

THEOREM 3 (Limit theorem). We have lim r .^ 00 = 

Tm(m— 1) Tm(m— 2) Mm (m— 2r+l) hm(m—2r) (/^m 1 hmo) 


Mm m /^m(m-l) Mm (m-r+1) Mm (m—r) Mm([m/2]+l) Mm [m/2] 

(Mm (2r—m) Mm (2r— m— 1)) (Mm m Mm (m-1)) 


Mm (m-r+1) Mm (m—r) 


Mm 1 Mm 0 


= $ (m) , (3.50) 


so that 


lim lx 2 - — 


(m,m) 


■ X + 1 X...X X + 


2 , 


r (m,l) 


X + 1 


= lim x 2 — T 


\ ( Af 

’’r -I I 2 . ‘"r , i 1 2771 






■ X + 1 X...X X + 




x + 1 = x zm + ... + x + l, 


T r 

where we note that as - fj, mk ) = (cj) mj - 2 - (0 mfc - 2)) = - <j> mk ), we 

can replace fi m i by f) m i in the above expression for \E r [ T "' ) . 

COROLLARY 1 (Convergence corollary). For 1 < j < rri, let a m j = Umi/hmj, 
so that d m i = 1, cr m 2 = (0 m l + 1) 1 ; and for j > 2, 0 < |(T m j| < 1. Let = 
4>mjt — 4>m(j-i)t> so that a[ m ’^ is the coefficient of [iff t in the closed form expression 


for given in S3.391 ) and f 3.48f >. Now define 

4m (I (f>mj — 0m. (j— 1) | T |</>m k 0771 (k— 1) I) — 2|0m j ~ 0 m (j — 1) | |0m k ~ 0m (k— 1) I 

|0mfc 0m (fc—1)| 2 

and let r k be the least positive integer satisfying 


1 


077 1 T 2 
Then for r > r k we have 

& 


rk 


< 


|0mfc 0m(fc—1)1 

2 (4m — |0mfc 0m(fc— 1)1) 


iT-f m,k) (m,k) 

V a) 


— j k | &m 2 1 j 
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so that approximates a'( n ’ J> /cL\ n ' K> , with remainder < 2Bj k |o‘ m , 2 | r . 

Furthermore, let r’ be the maximum of all the rk, B' be the maximum of all the Bjk, 
and for e > 0, let r* be the least positive integer with r* > r' satisfying 


( m,j) / ( m,k ) 


2B'y/m \(Tm2\ r <6. 

Then I'd/)- 7 ” — $( m )| < e for r > r*, and (by the standard definition in Section 1) the 
vector sequence is a sequence of vector convergents to the limit point 

When j = m — 2u + 1 with k = m — u + 1 = j + u, as described in | ITJfAty, or 
similarly as in \3.fO() , then |a^ m ’' 7 = \4> m u\, and the above simplifies to 

Am (\(f>mu\ + 1) OIJ , 

Bjk = —- 1 - \ — 2\<Pn 


so that for r > r^ we have 


?( m J) 


Am,j) 


a (m,k) 


< 


\tymk 4 > m(k— 1)1 


8 m ( \(fmu\ + 1) 
\ ( t > mk 1)1 


- 2|<^>„ 


"ml 


+ 2 


COROLLARY 2 (Continued fraction corollary). For each sequence of convergents 
to the limit point a!( n '^ /a( n ' k ^ generated by /^ m,k \ there will exist a corre¬ 

sponding sequence of simple continued fractions, which will converge to the simple 
continued fraction expansion of the limit point itself. 

Proof of Lemma 3.2. From (13.231) and (13.241) in the proof of Theorem 1, we can write 


Fmt PjiFmt) — l^mt ( Uj ( COS 


2irt 


2m + 1 

2 / irt. \ 


+ Uj -1 ( cos 


2irt 


2m + 1 


and using the identity \x m t = —4sin 2 J, yields 

irt 


HmtPjiUmt) = -4 sin 


2m + 1 


sin (27 r(j + l)t/(2m + 1)) sin (2irjt/(2m + 1)) 


sin (2irt/(2m + 1)) 


sin ( 2irt/(2m + 1)) 


cos (irt/(2m + 1)) 
-l 


sin 


_ — 2 sin (irt/(2m + 1)) ( . f 2irt(j + 1) 

\ 2m + 1 

nt(2j + 3) 


+ sin 


= cos 


irt 


2m + 1 


cos 


= cos 


irt 


2m + 1 


T' 2 j+:i I cos 


2m + 1 
irt 

2m + 1 


— cos 


2t Ttj 
2m + 1 

nt(‘2j ~ 1 ) 

2m. + 1 


- T 2j -i ( cos 


7 it 


2m + 1 


(3.51) 

by the definition of T n (x). Applying the relation T n+ \(x) = 2 xT n (x) — T n -\(x) from 
(11.61) of Lemma 1.1 and cancelling, we obtain 


= 2 cos 


T 2j+3 
irt 


cos 


irt. 


2m + 1 


- T 2 j- 1 ( cos 


irt 


2m + 1 


2m + 1 


T, 


2j+2 


COS 


irt 


2m + 1 


- T 2 j I cos 


irt 


2m + 1 


(3.52) 
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and substituting into (|3.51l) then gives /j, m t Pj(hmt) = 

T2j+2 {cos ^2^7+1") ) ~ T ^ ( COS ( 2m + 1) ) = ^ m U+pt~Tmjt = ( l ) m(j+i)t— ( l>mjt, 

which is (13.4411 . 

For (I3.45p . we write 

(fomrt (jfim (m—r+l)t (m-r)t) Pr (4*1711 /2)(T m _ r+1 (0 m t/ 2) Tm—r (0m t/ 2)) , 

and using the identity 2 T m {x)T n (x) = T m+n (x)+Tj m _ n | (®) from (11.611 of Lemma 1.1, 
after cancellation then gives us 

2 (T m _2r+1 (4>mt/ 2) — T m -2r (<f>mt/2) + Tm+1 (0 m j/2) — T m 2)) 

= 4 > m (m—2r+l)t 4 > m(jn—2r)t 4“ fmmt = 4 > m(m—2r+l)t 4 > m(m—2r)ti 

as required. Similarly we deduce (|3.46l) . 

□ 


Proof of Theorem 3. From (13.391) of the Corollary to Theorem 2, we have 


&(m,j+ 1) 

O' rj-> 


m m j 

'y ^ (fj’mt) Pj ihmt ) = y ] | (l^mt hj k) i 

t =1 t= 1 k =1 


and we consider 

^-(m,m-2fc+l) 2fc+l) / \l-r p / \ 

_ ^’r _ 2-/t= 1 1 Tint) *m—2k \Hmt) 

jrfam-k+l) - ^(m,m-2fc+l) ~ ^ ( /imt ) 1 " r P m _ fc (/w) ’ 


by (|3.44l) of Lemma 3.2. 

In the numerator sum above, for large positive values of r, the (/U m i) 1_r factor 
will dominate as the fj, m t are ordered in terms of increasing absolute value. Hence, as 
r —>• oo the above expression will converge to the ratio of the coefficients of (n m i) 1-r 
in the numerator and denominator. Therefore we can write 


lim 

r—>oo 


-(m,m—2/c+l) ' 




lim 

r—>oo 




(m,m—2/c+l)' 




(m,m—k+ 1) 


lim 

r—> oo 


| (/^ml) Pm—2k{h'ml) 

\ (/bnl) Pm—k 


_ (hml)Pm—2k{hml) _ 4 > m(m—2k+l) 4 > m(m—2k) _ ^ 

{h"m l) P-rn—k ( t > m(m—k+ 1) — fim^m—k) 

by (13.4411 and (|3.45j) of Lemma 3.2, and hence the result. The polynomial identity 
in the limit then follows immediately . 

Regarding the first Corollary, let ( , so that by f)3.48[l 




_ J m d) 

— u, 


+ a 


(m,j) 


(7, 


m 2 


+ ... + a 


(mj) 

m 


a 


r 

mm 



a l 



A m J) 

a l 

gr(™,k) 
c_x f 

(m.k) 

a\ 


jgy(m,k) 

(m.k) 

«i 


and 
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— a 


1 


(m,k) 


V 1 


- 


(m,k ) ( m,k)\ (m,j) 


— a 


a 


( m,k ) / (m,fe) 


1 


1 


+ 


(m,k) 

jTr - a 


Now for 1 < j < m, we have 2 




< 4m — la 


( m >i) i 


i 


and 


jrr (mj)_ a (mj)l < ( |{ 


(mj) I 


+ . . . a, 


(™.ii 


cr: 


m 2 b 


so that 


^( m d) 


,(mj) 


< 


4 m ’ fc) |(|af J) | 


H - ... Cl' 


eyoj- 

( m j) 


(m,k) 


(m,k) 


+ |aS mJ) 


I(i4”‘ , ‘ > i 


+ 


(m,k) | 


(m,k) | /1 ( m,k ) i 


a. 


(m,k) | 


i l (wijfc) 
+ . . . Ofn * 


<T, 


r—1 1 

m2 I 


<T, 


^r-li 
m2 I 


Hence for r > with ras defined, the results follow. 

Regarding the vector distance |V&r — by the definition of B' and r' it 

follows that for any specific values for j and k with r > r' , the remainder term will be 
< Taking the square root of m times the square of this bound to obtain 

the Euclidean distance, we see that the remainder term is now < B'^/m\a m 2 \ r ~ l • As 
\(7 m 2I < 1, we can always find such an r* > r' satisfying the inequality, as required. 

To see the second Corollary, for an irrational number a, the simple continued 
fraction algorithm begins with the largest integer not exceeding a and then proceeds 
to generate what can be thought of as one infinite branch of number plus fraction. 
For the (m + 1) term recurrence relation that generates the sequence of numerators 
and denominators in /,^ ; ! r rn ' k> . at each step the continued fraction branch will 

have number plus m — 1 fractions, so at each step the number of branches will 
increase by a factor of m — 1. For m > 3, it is expected that in general the resulting 
sequence of rational numbers will converge at a slower rate than that for the simple 
continued fraction expansion of a( n '^ / a( n ’ k \ as guaranteed by Dirichlets theorem 
(11.131) . However, the resulting sequence can itself always be translated into a sequence 
of simple continued fraction expansions, and by the theory, this must converge to 
the unique simple continued fraction expansion of the limit point ag 

required. □ 

Definition (of unlaced Fibonacci sequences). For 1 < j < m, denote by 

the rth term in the jth unlaced Fibonacci sequence of dimension m, defined in terms 

of the interlaced Fibonacci sequences of dimension m such that 




(m,j) 


mq-\-2r—l 


mq-\-j ’ 


1 < r < 


m + 1 


and 


^ mq-\-2r ^ mq+j 


1 < r < 


m 

lTj 


so that the unlaced sequences of Fibonacci numbers correspond to the ordered union 
of the ordered sets 

OO 

j* _ I I ( G?(m, 1) 1) «-(m,2) 2) ®-(m,[m/2]+l) \ 

^ r J n i^n i^n i^n i • • • " n J • 

q =o 

n=qm+j 
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Here the ordering mimics that displayed in the construction of the convergent vectors 
terms where as before, we define the unlaced numerator terms to be 

the non-reduced numerators of the terms 


LEMMA 3.3. For 1 < o, b < m, a ^ b, and with the above construction, the ratios 
of the terms of the jth unlaced Fibonacci sequences of dimension m, given by 


* 

o' , 


mt+a 




(™,j) 


mt+b ) t =0 


n Ar( m d) 
^ r mt+a 

„ jyirnj) 
tyr mt-\-b 


t =0 


are a rational sequence of convergents to some real limit points, including those of 
the form 2 cos • 


Proof. The proof follows immediately from the construction of the unlaced Fibonacci 
sequences of dimension m from the interlaced sequences of dimension m. □ 


Example (Three-dimensional interlacing case). By construction, the interlacing 
Fibonacci sequences of dimension 3 are given by 


/ Jf’ 1} 




(3,2) 
r(3,3) 




(3.1) 

(3.2) 

(3.3) 


^(3,D \ 

jr(3,2) 


r(3,3) 


y e ^2 >,OJ 

and thereafter for j G {1, 2, 3} by 



-3© 

-I© 

-id) 





r+3 


= -2& 


r+2 




*(3j)_ *(3 ,j) 

r+1 r > 


so that 



(3.53) 


/ *( 3 4) 

/ ^r+3 

*(3,1) 
tS r+ 2 

*(3,1) 

^r+1 

\ 


/ *(3,1) 

^r+2 

*■(3,1) 

^r+1 

*(3,1) 

ex 7* 

\ 

/ -2 

1 

0 \ 

jr( 3 > 2 ) 

^r+3 

*(3,2) 

^r+2 

*(3,2) 
^ r+1 


= 

*(3,2) 

^r+2 

*■(3,2) 

^ r+1 

jr(3,2) 



0 


\ *( 3 ,3) 

\ ^r+3 

*■(3,3) 

^r+2 

*(3,3) 

'A’+l 

j 


*(3,3) 

\ "A-+2 

*(3,3) 

^r+1 

*(3,3) 

tA t* 

J 

\ 7 

0 

0 ) 


By (13.331) of Theorem 2, the generating functions for the sequences , 7?i 2 ' ) and 
&r^ are given by 


*( 3 ,i) . 7(x 2 + 4x + 3) 

x 3 + 7x 2 + 14.x + 7 


3—2x+2x 2 


17 x 3 
——h 


22x 4 

7 


29x 5 | 269x 6 357x 7 | 

7 ' 49 49~ + 


jr(3,2) 

*■(3,3) 

o' rp 


7(x + 2) 0 0 9 37x 3 49x 4 

: ——Gt- -L - = 2—3x+4x 2 --—- 

x 3 + 7x 2 + 14.x + 7 7 7 

7 , „ „ 9 29x 3 39.x 4 

x 3 + 7x 2 + 14x + 7 7 7 


65x 5 , 604x 6 802x 7 , 

7 ' 49 49 

52x 5 | 484x 6 643x 7 ( 

7 ' 49 49~ + 


and by (13.341) JFr 3 ' 1 ' 1 and can be expressed in terms of j^ 3,3 ^ as 


jK 3 - 1 ) = 3jH3,3) + 4^ r (3 ’ 3) + JF 


(3,3) 
r —2 ’ 


J?( 2 ) — 2^"( 3,3 ( + f 2‘( 3,3 ) 


r—1 


The first few terms of the three unlaced sequences *,!pJ: rn ’ J ' > . with j G {1,2,3,}, 
and their corresponding non-reduced numerator integer sequences *.jV} 3 ' 3 ' 1 are given 
below. 
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Hence, from the definition of and by Theorem 3, we have 


* 

lim r = lim 


^■(3,2) * 27T 

——— = lim - 7 ^—— = 2 cos — = 1.24698 ... 


r( 3 >3) 


r—> oo * 




(3,1) 


.3)— 1 


r—>oo * <gj-(3,2) 
^3r-l 


i—>oo * 




(3,3) 


3)—1 


* ®"(3>1) 

r ^3r-2 
hm 

r—» oo 




(3,1) 


*+7 


= lim ■ 

r—» oo 


3r 

* ®"( 3 d) 

lim ^ 

r—> oo 


(3,1) 


*+7 


(3,2) 


3i—2 




(3,2) 


3r 


= lim ■ 

r—>• OO : 




3i—2 




* ^j-(3,3) 
- 1- ^3r-2 

i—»oo *^-(3,3) 


3r 


(3,2) 


3r—1 


(3,2) 


= lim 


707 


(3,3) 


3i—1 


3r—2 




3i—2 


—2 cos — = 0.445042 . 
7 


6n 

-2 cos — = 1.80194. 
7 


which is similar to the ratios of consecutive Lucas numbers or Fibonacci numbers 
converging to —2 cos or their reciprocals to 2 cos nf. We note that due to these 
sequence structures, the rational sequences can be replaced with the non-reduced 
numerator integer sequences without affecting the ratios. 


Example (Of cyclotomic approximation when m = 5), 
Euclidean distance regarding the 20th convergent when m 


In consideration of the 
= 5, we have 


^20 ^5 




(5,4) ^(5,2) 


20 


20 




(5,1) 


20 




(5,3) 


20 




(5,5)' 


20 




(5,5) ’ ^(5,4) 


20 


20 


& 


(5,3) 


20 




(5,2) : 


20 




(5,1) 


— (4>5l,<t>5 2,4 , 5 3, ( t , 5 4:,4 > 5 5) 


20 


( 42951850444254470 35685687021511133 4434370056070408 46738310388496383 25528481467235249 \ 
V25528481467235249 ’42951850444254470 ’ 15579436796165461’ 35685687021511133’ 13303110168211224 / 




, 2 cos 




2 cos 



< 1(T 13 . 


In terms of polynomials, this yields (to two decimal places) the remainder term 



*(5,4) 

■+0 

*( 5 , 5 ) 

•^20 





*( 5 , 1 ) 

£_ 20 _ 

*( 5 , 3 ) 

•^20 




*( 5 , 3 ) 

■^20 

*( 5 , 2 ) 

•^20 




*( 5 , 5 ) 

■+0 

*(5,1) 

•^20 



— (.T 10 + x 9 + x 8 + x 7 + (E 6 + x 5 + X 4 + x 3 + X 2 + X + l) 

= -1.48 x HT 13 x + 1.10x 10“ 2 V-2.96x 10 _13 x 3 + 2.20 x 10 _2 V - 2.96 x 10 _13 x 5 
+2.20 x 10 _2 °x 6 - 2.96 x 10 _13 x 7 + 1.10 x 10 _2 °x 8 - 1.48 x 10 _13 a: 9 . 
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Definition (Of Fleck’s and Weisman’s congruences). Let pbea prime and a be an 
integer. In 1913 A. Fleck discovered that 

(—l) fc (N') = 0 ("mod p\- p _1 -I , (3.54) 

k=a (mod p) ^ 


for all positive integers N > 0. In 1977 C. S. Weisman [23] extended Fleck’s congru¬ 
ence to obtain 


E (-+ 

k=a (mod p a ) 


= 0 (mod p w ), 


co = 


N—p° 


0(p“) 


(3.55) 


where a, N are positive integers > 0, N > p “ , and <f> denotes the Euler totient 

function. When a = 1 it is clear that (13.551) reduces to (I3.54D . 

We define the Fleck numbers , ${N,a (mod n)), such that 


$(N, a (mod n)) = ^ (3.56) 

k=a (mod n) 

These sums have many well known properties m, pe] such as 

^ / j\T\ 

nd{N,a (mod n)) = ^(-l) fc ( J £ j k ~ a = ^ 7 ~ a (l - 7 )^, (3.57) 

k= 0 ^ ' 7"=1 7 n =l 

from which we can deduce the recurrence relation 


$(N +l,o (mod n)) = ${N, a (mod n)) — ${N, { a — 1) (mod n)). (3.58) 

By modularity we also have ${N,a (mod n)) = $(N, (a + n) (mod n)). 

In the following theorem we give new expressions for the renumbered Fleck num¬ 
bers in terms of our polynomial functions. 


THEOREM 4. Let r be a non-negative integer and m be a positive integer, so that 
n = 2m + 1, is odd. Then the numbers in the sequences and are given 

by the alternating binomial sums 



= n$(2r + l, {r + j) (mod n)), 


(3.59) 


A m d) 


= n 




2r+ 2 


\r + j + 1 + an 


nff{2r + 2, (r + j + 1) (mod n)), 


(3.60) 

and so when n is a prime power, they satisfy an analogue of Weisman’s Congru¬ 
ence fEEf. 

The renumbered Fleck numbers can be written as n${2r + 1, (r + j) (mod n)) = 


^ ( {Amt) -fj — 1 {Amt) — ^ ^ {Amt) Vj—1 ( ‘fmt ) ^ ' {Amt.) {tymjt f’m (j—l)t) • 

t= 1 t= 1 t= 1 

(3.61) 

Similar expressions exist for the sequence terms , corresponding to the renum¬ 

bered Fleck numbers ${2r + 2, {r + j + 1) (mod n)), formed from the even rows of 
Pascal’s triangle. 
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Proof of Theorem f. These relations are obtained by Lemma 4.2 of m and the 
closed form expressions (13.4811 . and (13.4911 of Lemma 3.2. □ 


Example. When j = m = 5, so n = 11, and r = 8, we have 

11 E(-u‘ 




17 

2 + 11a 


5 5 

= (l“5 1 ) 9 Pi (fist) = J2 ^5t) 8 0^55 1 ~ 054 1 ) = 11^(17, 2 (mod 11)) = -11 2 .204. 

4=1 4=1 

Remark (to Theorem 4). Using the recurrence Q m (x) and its inverse recurrence, 
we can construct a similar family of sequences, which at negative indices correspond 
to (2m) times the Fleck numbers obtained using the even modulus (2m). 


There are many geometric interpretations of the Fibonacci numbers and the 
Golden Ratio, including the triangles inscribed in a circle representation, given by 
J. Rigby in m- In the following theorem we establish geometric relations for our 
rational interlacing Fibonacci sequences (and % ), with the ratios of di¬ 

agonal lengths between the vertices of the regular n -gon inscribed in the unit circle. 


THEOREM 5. Let n a positive odd integer, with n = 2m + 1 , and v n o, ■ ■ ■ v n ( n -i) 
the n vertices of the regidar n-gon H n , inscribed in the unit circle. Let d n k be the 
signed distance from vertex v n o to vertex v n k, so that d n \ = 2 sin ( 7 r/n) is the side 
length of H n , d n k = 2sin ((nk)/n) is the signed length of the kth diagonal of H n , 
and d nn = 0. Then working (mod n) the second subscript of d n k, we have 


= (-i) r - x x. = w 1 E 

4=1 4=1 (2sin^ nt 


JJ-(rnJ) _ (_-Qr—1 lS)t 

t 1 

and 

m 

= (-ir‘E 


d. ... 2sin 


(2j-l)7T4 

2m+l 


2m+l 


2r—1 ’ 


(3.62) 


_X d n d n (2j— 1 )4 


4=1 


j2r—1 
l nt 


(-1 r‘E 


2cos[Jgt. 


4=1 (2 sin 


7r4 

2m+l 


2r—2 ' 


(3.63) 


Hence for r a negative integer, these sums over the signed diagonal lengths correspond 
to integers which are the renumbered Fleck numbers multiplied by n. 


Proof of Theorem 5. Combining the definitions of the diagonal distances and ratios 
with the closed form trigonometric expressions for given in the Corollary to 

Theorem 2, we obtain (13.621) . 

Applying the trigonometric identity sin x —sin y = 2 cos P-Fll sin — to ^4 m d+ 1 )_ 
&r\ we find that {d n ( 2 j+i)t—d n ( 2 j-i)t)/dnt = 4>mjt , and hence (13.631) Vr E Z. □ 

Remark (To Theorem 5). Defining the diagonal to side length ratio r n k, such 
that r n k = d n k/dni, we note that the sums over ratios of signed diagonal lengths, 
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corresponding to ratios in the unit n- gon, are closely allied to Steinbach’s unsigned 
diagonal and ratio product formulae (see m and 0 ) given by 

e-i £-1 

\dnkdn£\ = \d n l| \d n £-(- 2 j+l) | ; \ r nk r n£\ = ^ ^ I r n (k— £+27+1) I; 

j=0 j=0 

where we take (A; — •£ + 2j + 1) (mod n). 

In the following theorem we examine the non-reduced numerator integer inter¬ 
lacing Fibonacci sequences of dimension rn. , in more detail. 


THEOREM 6 (Integer sequence theorem). Letp\,p 2 , ■ ■ -Pt be all the distinct prime 
factors of n = 2m + 1, so that for each pt we can write pi = 2q l + 1. Then for 
1 < j < m, and r > 0, we have 





g Z- 


For with r < 0, and n having at least two distinct prime factors, let 

; an d when n = p a is a prime power, let 

= ) ■& = -a- 

Then with n = jq 1 ... pf*, we have for r = 0, —1, —2, —3,..., that the sequence terms 
, are t/ie renumbered Fleck numbers, and when n = p a is a prime power, the 
sequence terms are the (renumbered) Fleck quotients discussed in ll8f . 


—2 r + 1 — p° 


pa _ p. 


a— 1 


COROLLARY. W7ien n = p, a prime number, we have 

{pL^Jjr("M')\ +0 ° = {^Ki)j + °° ) so that eZ, V r eZ. 

I J r=—oo l J r=—oo 


Proof of Theorem 6. The case n = p, a prime number, was proven in Lemma 7.5 
of flUll . An expansion of this argument then leads to the deduction that for n = 
P^pT ■ ■ ■Pt *! an d r positive integer values, the power of pi in the denominator 
of the sequence terms never exceeds |_{r — l)/qi\. Hence we can treat each prime 
number individually, replacing the denominator m in the exponent of pl4 r ~ 1 )/ m J, 
with q t for each prime factor p t . It follows that multiplying together the product of 
the prime numbers to their individual exponents yields an integer sequence. 

By Theorem 4, we know that the sequence terms are the renumbered Fleck 

numbers, and so are already integer values and satisfy an analogue of Weisman’s 
congruence. The property relating to Fleck quotients when n = (2m + 1) = p. a 
prime number, then follows from the exponent so that 



1 

to 

-i 

1 

to 


1 

to 

*i 

1 

to 


—r — 1 

a=l 

p — l 


2 m 


m 


as required. 


□ 
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Example. We give the table for the integer interlacing Fibonacci sequence 
when m = 3, and r € [—6, 7]. 


j/r 

-6 

-5 

-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 

6 

7 

1 

-35 

66 

-18 

5 

-10 

3 

-1 

3 

-2 

2 

-17 

22 

-29 

269 

2 

26 

-47 

12 

-3 

5 

-1 

0 

2 

-3 

4 

-37 

49 

-65 

604 

3 

-13 

22 

-5 

1 

-1 

0 

0 

1 

-2 

3 

-29 

39 

-52 

484 


Remark (to Theorem 6). When n = 2m+l has repeated prime factors, it appears to 
be the case that particular sequences out of the m sequences have further divisibility 
properties. One area for future investigation would be to derive a concise formula 
that describes these patterns. 


LEMMA 3.4 (Sum of diagonal products lemma). Let H n be the regular odd-sided 
unit n-gon H n , as defined in Theorem 5, so that n = 2 m + 1 > 0, H n has signed 
diagonal distances given by d n k = 2sin((7r k)/n). Then we have 


m 

^ ^ + (2j—l)u dn (2j—l)v 
3 =1 


0 if ufi v 

2m + 1 if u = v 


Proof. Applying trigonometric identities we can write the sum as a difference of sine 
ratios, such that 


^2 sin 
3 = 1 


iru(2j - 1) 
2m + 1 


2 sin 


( 2mir(u—v) 
1+2 m 

2rn + 1 / ( rc{u-v) 

6,111 1 l+2m 


sm 


sm 


( 2rmr(u-\-v) \ 
1+2 m ) 

7 t(u+v) \ ’ 

l+2m I 


thus eliminating the j variable. Using sin (7r — x) = — sin a;, we see that the two 
terms cancel out if u v, and when u = v, we deduce the result via the identity 


sin 2 (£ix) = 1/4(1 + 2 N — cscxsin[x(l + 21V)]). 


□ 


THEOREM 7. Let n = 2m + 1 be a positive integer. Then the sequence terms 
and&r" 1 '^, obey the relations 




(m, 1) _ 1 


2 r 


n 


E * 


r 


ST. 


(m, 1) _ 1 


2r+l 


3 =1 


2 n 


m, 0 ) 


+ 2 E ( Cg{rn,j) 

3 = 1 


so that for all r € Z we have 


-1 


m / m 


^2 ++ 2, +1 = — r i^rnjt - fim(j-l)t) 


t= 1 


71 j=l \t= 1 


and YT=i (hint 


,—2r 


Y l m / m \ ■ 

= I 4 (/^mt) + 2 ^ ^ ^ (hint) jt "h fim (j—l)t) y 

Hence each of the sequence terms and + mj ' 1 can 6e written as an integer 

linear combination of sums of square integers. 
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COROLLARY. Let n = 2m + 1 > 0. Then we have the Fleck number relations 

m 

5(4r + 1, (2r + 1) (mod n )) = — ^ (5(2r + 1, ( r + j ) (mod n))) 2 , 

i=i 

25(4r + 3, (2r + 2) (mod n)) = 5(4 r + 4, (2 r + 2) (mod n)) 

m 

= (5(2 r + 2, (r + 1) (mod n))) 2 + 2 ^ (5(2r + 2, (r + j + 1) (mod n))) 2 . 

3 = 1 

T/ie corresponding combinatorial identities (not given in fffj) are 

\ 2 


m / oo 


E (-^TAiTJ-E E (-i r +i+ ‘LlT + ln) 

i =—oo j=l \a =—oo , 


2 E = E 

2 


a=—oo 
m / oo 


E (-i)’' + 1 + “C’,?J + 2 E E <- 1 ) rW+ 1 + “U|«+a„) • 


\a =—oo / j =1 \a=—o o 

We note that the above identities appear to be true when n = 2m. 


Proof. Writing the first display of the Theorem in terms of the signed diagonals of 
the unit n-gon H n and rearranging, gives us 

m m 

^^(dnn) 4r+2 {d n (2j-l)u) 
u=lj =1 

mm m 

+2^ ^ (d nu ) 2r+1 (4ti) 2r+1 ^ d n (2j-i)u d n (2j—i)v 

u= 1 v=u +1 

By Lemma 3.4, for u ^ v the inner sum is zero, and by (13.481) the result follows. The 
second display can be obtained similarly. 

To see that each sequence term can be written as an integer linear combination 
of sums of square integers, by the preceding results we can write each term # 2 ™' *^ 
as a sum of m squares, and each term ^ 2 r+i as a sum °f 171 squares plus twice a 
square. By (13.341) and ()3.35l) . each of the sequence terms and , can be 

written as an integer linear combination of the , and hence the result. □ 

Example (Of sum of squares representation). In consideration of the non-reduced 
numerator A^g 6 ’ 3 ^, of the sequence term J^g’ 3 ^, by (I3.34|) can write 

-^ 6,3) = 29226191 = - 5A^ (6,1) - ^ (6,1) 

= 5 x 7480420 - 5 x 1713705 + 392616 
= 5 x (1505 2 + 1421 2 + 1245 2 + 1010 2 + 702 2 + 365 2 ) 

-5 x (702 2 + 645 2 + 543 2 + 411 2 + (2 x 365 2 ) + 260 2 + 84 2 ) 

+344 2 + 327 2 + 283 2 + 234 2 + 159 2 + 85 2 , 

so that after simplification, this method enables oT)g ’ ' to be written as an integer 
linear combination of 16 square integers. 
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4 Properties of the Polynomial Generating Functions 

THEOREM 8. With Pq{x) = 1, P\(x) = 1 + x/3, Qi(x) = 1 + x/2 and Q 2 (x) = 
x 2 /4 + x + 1/2 ; the polynomials P m (x) and Q m (x) respectively satisfy the three-term 
recurrences 

Pm+i(x) = {x + 2 )P m (x) - P m -i(x), (4.1) 

Qm+i(x) = {x + 2 )Q m (x) - Qm-i(x), (4.2) 

the ordinary differential equations 


x(4 + x)P'f(x) + 2(x + 3 )P' m {x) - m(m + 1 )P m (x) = 0, 

x(4 + x)Q'/ n (x) + (2 + x)Q' m (x) + m 2 Q m (x ) = 0, 
and for integers k,£, the explicit orthogonality condition 




4/2 


l/^mx) {4+x)1/2 


dx = 2m8tki £,k > 0 


r-O 


Qe(x)Q k (x) 


y_ 4 ,x i / 2 (4+ x) 1 / 2 
where bn k is the Kronecker symbol. 

Let 

f° P m (x)x s ~ 3 / 4 


dx = —2TTidnk, k f 0, 


m p = r p n Ax}.i"' - 

m ( s ) ~ L A (4 + x) 3 / 4 


dx, Re s > —1/4, 


r Qm(x)x S ~^ , 

MS(S) = /_4 (4 + x)3/4 

TTien up to normalization, these Mellin transforms have the form 

r ( s + i) 


M£(s) = (—l) s+5 ^ 4 4 s 4 _m_1 r(l/ 4)p m (s) 


M%(s) = (—l) s+3 / 4 4 s_1 r(5/ 4)q m (s) 


r( s + ^i)’ 

r(s-i) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
(4.10) 


T(.s + m) 

COROLLARY. Closed form expressions for P m (x — 2) and Q m (x — 2) are given 
by P m (x - 2) = 

Vx + 2 \ r r-= -\ m , A , V® + 2 ' 


■>—m—1 


(( i _^) ( ._^ r+ ( 1+ 


and 

Qm(x- 2) = 2~ r 

where for r < [m/2] we also have 


Vx^2 

-V / ^) m +(x + V^) m ) 


(4.11) 


(4.12) 


P m {x) = £(-l ) r ( r ) (x + 2) r ~^P m - r -j(x), (4.13) 

3 =0 
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(4.14) 

(4.15) 

(4.16) 


and 

Qm(x) = £(-l y( r \x + 2) r - j Q m - r - j (x). 

3=0 

We have the orthogonal polynomial relations 
[0 

/ X 1 ' 2 (4 + x)' - 1/2 a ( x)(any polynomial of degree < k)dx = 0. 

[0 

/ x~ -V 2 (4 + x r ,,2 Q t (x)(any polynomial of degree < k)dx = 0. 

The polynomial factors of M^(s) and M„ r (s) satisfy the functional equations 

Pn(s) = ±p n ( 1 - s), q n (s) = ±q n ( 1 - s) (4.17) 

and have zeros only on the critical line Re s = 1/2. 

Proof of Theorem 8. The three term recurrences for P rn (x) and Q m (x) in (14.11) and 
(14.21) follow from the Legendre function expression for P m (x) given in (13.221) . and 
using ID [p- 99 or 247 or 295] with a = (3 = —1/2 the Jacobi polynomial Pn°^\x) 
relation 

_ (2m)(m-l)! ( _ 1/2 ,_ 1/2) / x\ 

X) ~ (l/2) m m V i+ 2/' 

We sketch the details for Q m (x). A recurrence for the polynomial Pf ’ ' (z) is 

(n + l)n(2n - l)P^/ 2 ~ 1 ^ 2 \z) 

= n(2n - l)(2n + 1 )zPj~ 1/2 '- 1/2) (z) - {n - l/2) 2 (2n + l)pjfl\ /2 '~ 1/2) (z), 

and using the change of variable z = 1 + x/2 and (2.5.14) of [l], applied to the 
polynomials Q m (x), we obtain 

Considering now the ordinary differential equation satisfied by P„(z), 


. n, d 2 u du 

(1 -*V- 2 7fe + 


v{v + 1) — 




u = 0, 


1 — z 2 

and an elementary application of the chain rule, we find 

x(4 + x)P^(x) + 2(x + 3 )P' rn (x) — m(m + 1 )P m (x) = 0, 
and hence (14.3|) . 

By using an integrating factor x 3 / 2 (4 + x) 1 / 2 , the differential equation for P m (x) 
may be written as 

£ 3 / 2 (4+£’) 1 / 2 P" 1 (x)+2(x+3)£ 1 / 2 (4+x)~ 1 / 2 .P) n (x)—m(m+l)x 1,/2 (4+x’) _1 / 2 P m (x) = 0. 


We then obtain 
d 

dx 


[a?/ 2 (4 + x)V 2 P^(x)] = rri(rn + 1) x 1 / 2 (4 + x)- 1 / 2 P m (s). 
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Writing this equation for Pfc(x), multiplying the P m {x ) equation by Pfc(x), and the 
Pk(x) equation by P m (x) and subtracting there follows 

P*M A[x 3/2( 4 + x) 1 / 2 .P m (x)] - P m (x)-A{x 3/2 (4 + x^fx)] 

= [m(m + 1) - k(k + l)]x 1,/2 (4 + x)~ 1 // 2 ]Pfc(x)P m (x). 

Thus 

{x 3/2 (4 + x) 1 / 2 [P fc (x)P/ n (x) - P m (x)P^(x)]} 

= [m(m + 1) — + 1)]x 1//2 (4 + x)~ 1 // 2 ]Pfc(x)P m (x). 

By integrating between x = —4 and 0, we obtain the stated result (14. 5 j) . 

The orthogonality of the sequence {P m (x)} m >o allows the development of integral 
transforms with zeros only along vertical lines in the complex plane and hence (14.71) 
and (14.91) . 

For (14.4p we use the hyperbolic trigonometric function analogue for Chebyshev 
polynomials to obtain 


Qm(x') 


2 cosh 


^ 2 m sinh 1 



22F1 



The differential equation for the Gauss hypergeometric function 2 F 1 (a, 6 ; c; z), 


z(l 


yd 2 ydy 
z) lgI + lc-( a + b + l)z] Tz 


aby = 0, 


becomes for Q m (x) 


d 2 y (2 + x) dy m 2 _ 

, 9 H- t~. --1-- ry — 0 

ax z x(4 + x) dx x(4 + x) 


and the result follows. 

In consequence, the family {(J m (x)} m > 1 is orthogonal, and with the integrating 
factor yfxyj A + x, the differential equation may be written as 

,2 


d_ 

dx 


( vwiT *!0 


m 


y/xy/A + X^ 


The integrating factor is obtained as the exponential of 

f (2 + x) , 1 . r , A 

J x(4 + x) 2 L v 

We then obtain the orthogonality relation (14.61) [the steps being omitted] 

Qm{x)Qk(x) 


,1: 


;dx = —2iri5 m k, k / 0. 


/_4 x x / 2 (4 + x ) 1 / 2 
Accordingly, we have a (generalised) Mellin transform 

r ° Qm(x)x S ~ 5 / 4 


M%{s) = J_ 


-4 (4 + x ) 3 / 4 


dx, 
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of the form (|4.10|) . so that 

M%{s) = (~l) 5 + 3 / 4 4 s ~ir(5/4)g m (s) F ( S ~ . 

i (s + m) 

In the Corollary, (14.111) and (14.121) are obtained by solving the recurrences in (14.11) 
w, whereas (14.131) and (14.141) arise from iteratively applying the recurrences. 

The last part of the Corollary follows from properties of orthogonal polynomials. 

The proof that the polynomial factors of M^(s) and M r „ (s) satisfy the functional 
equations p n (s) = ±p n (l — s), and q n (s) = ±<7 n (l — s), and have zeros only on the 
critical line Re s = 1/2, follows that given in |3j 

□ 


THEOREM 9. The polynomials P m (x ) obey the Christoffel-Darboux formula 


E Pm(y)Pm(x) 

m =0 


P n+ l{x)P n {y) - P n +l(y)Pn(x) 
x-y 


the confluent form of which (i.e., for y —» x) is 


(4.18) 


= P 'n+1 (x)Pn(x) - P n+l {x)P' n {x). (4.19) 

k =0 

They also satisfy the relation 

—= ,, , t {-(2m - 1 )P m -i(x) + [m(x + 2) - 1 ]P m (x)} , (4.20) 

ax x[A + x) 

and have the raising and lowering operators 

Rm = x(x + 4)-^- + m(x + 2) + x + 3, (4.21) 

dx 

L m = x(x T 4)-y- f- m(x T 2) 1, (4.22) 

dx 

such that R m P m (x ) = (2m + 3)P m+ i(a;), and L m P m (x) = (2m — l)P m _i(x). T/ie 
ordinary differential equation for P rn (x) can then be written in terms of these oper¬ 
ators. 

The polynomials Q m (x) obey the quadratic identity 


Qm( x ) — Qzm (x) + 2 , 


(4.23) 


have the generating function 


V { l^LQ m {x)r m = R-\l-r + R) 1 / 2 ^ + r + R) 1 / 2 , 

Z —' ml 

m =0 


where R = (1 — 2r — xr + r 2 ) 1 ^ 2 , and satisfy the differential relation 


. Qm(x ) 

dx 


2msin[mcos 1 (l + x/2)] 
\/—x(4 + x) 


(4.24) 


(4.25) 




















37 


Proof of Theorem 9. The polynomials P m (x) have the hypergeometric form 

2 


P„W = 2 F I (- m , m +l;|-^j = (2m+lv _ 


sinh 


/ / % 

(2m + 1) sinh -1 ' 


Hence their ODE may also be found from that of the 2 F 1 function. 
If we normalize such that Pp(x) —>• Pp(x')j\j2'xi^ so that 


/ 


0 

-4 


Pk(x) 


4/2 


-dx = 1, 


(4 + x) 1 / 2 

we obtain the Christoffel-Darboux formula of ()4. 18|) 

P n+ i(x)P n (y) - P n+ i(y)P n (x) 


22 Pm(y)Pm( 


X) = 


m =0 


X - 7/ 


and the confluent form of this result (14.191) 

n 

J^pf(x) = P' n+l (x)P n (x) - P n+1 (x)P' n (x), 

k=0 

then follows. 

When the relation 

(1 - z 2 ) 


dP 


- 1 / 2 , 


dz 


= m 


- P m-l( z ) - ™z p m 1/2 ( z ) 


is transformed to the polynomials P m (x), the result is 
dP m (x) 1 


dx x(4 + x) 


{-(2m - l)P m _i(x) + [m(x + 2) - l]P m (x)} 


which is (14.201) . The raising and lowering operators of (14.211) and (14.221) can then be 
deduced. 

From the application of linear and quadratic transformation of the 2 F 1 function 
we have the following. 


Qm(x) =2 1 + 


T? I 1 1 X 

2 F 1 -m, - - m; 


2 x + 4 


(4 + x) 
2 2m 


1 - 


4 + x 


2 m 


+ 1 + 


X 


4 + x 


2 m 


=21+ 


= 2 i i+ i 

As 


x \ 1/2 


_ . 1 lx 

2 ! i m, —im;-;- 

1 ’2 ’ 2’ x + 4 


1 x 


) 2 F ! I ^ + m, ^ - m;1 = 2cosh I 2msinh 1 I -y- 


^ , m m 1 x., , , 

g m (x) = 2 2 Fi ( ~ y , y ' 1 2 , ~^ 4 + x ^ ^ ’ 
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we find that Q m (x) = Q m / 2 [s(4 + x)]. As 


Q m (x) = 2(-l) m 2 Fi 1 + x/4'j , 

we determine that Q m (x) = (—l) m Q m (—4 — x). Furthermore, 
2^/n(2 m — 1 )! 


Qm{x) — 


(m — l)!r(m + 1/2) 

2^tt( 2 m — 1)! (x 


( x\ m ( 1 1 \ 

( 1 + l) 2Fl 1,-“' 2 - m; 1 ~ 2m; TT^I j 


(m — l)!r(m + 1/2) \4 


Fi I —m, - m; 1 — 2m; —- 


x 


With (a) n the Pochhammer symbol, we note the limit for j > 0 


lim 


(a + b)j 


a-^r—b (2 CL + 2b) j 2 ’ 

otherwise this ratio is 1 for j = 0. We then obtain a reduction of Clausen’s identity 
for the square of a special 2 Fi function, Q^(x) = Q 2m (x) + 2, which is (14.231) . 

To see (14.241) we identify Q m (x) in terms of Jacobi polynomials P^'^\x). We 
use ID [p- 99 or 247 or 295] with a = /3 = —1/2 and obtain 

O (x) ~ P(-V2-1/2) (l , 

J '(l/2) m m l i+ 2J’ 

which is the Gegenbauer polynomial case of C^°. 

A generating function for Jacobi polynomials is [I] (p. 298) 

OO 

F(z,r) = ^P^\z)r n = 2 a+f} R~ 1 {\ - r + R)~ a {\ + r + R)~P, 

n =0 

where R = (1 — 2zr + r 2 ) 1 / 2 . Correspondingly we find the generating function 


E 

m =0 


(l/2)r 


ml 


-Qm(x)r m = R~ 1 (l-r + R) l/2 { 1 + r + i?) 1/2 , 


as required, where now R = (1 — 2 r — xr + r 2 ) 1 / 2 . 

By P (p. 297) 

± p (- 1 / 2 -i/ 2 ) f n = n p (i/ 2 ,i/ 2 ), x = r(n + l/2)sin(ncos- 1 x) 
dx n 2 n_1 ^(n-l)!^!^ 


We then obtain 


d 2msin[mcos 1 (l + x/2)] 

=- 7 - 

dx yj— x(4 + x) 


which is (I4.25p . 

We note that in terms of Jacobi polynomials P m (x ) can be written as 


P - m! p( 1/2,-1/2) A + 

mW_ (3/2) m m V + 2/' 


□ 
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Remark (to Theorem 9 - Generalized raising operator and Rodrigues’ formula). It 
is possible to obtain a generalised Rodrigues’ formula for the polynomials P m (x), as 
we now present. Following the procedure of j2j we put R m = fi/fg 2 + h, where h is 
an arbitrary function and f\ and g 2 are functions to be determined. We find 


g 2 (x) = exp 


h(x) 
x(x + 4) 


dx 


x {m+ 3/2)/2 + (m+l/2)/2 , 


and 


x(x + 4 

h(x) = —= ex P 

92{x) 




h(x) 


x(x + 4) 


dx 


x ~ m / 2+1/44^-m/2+3/4^ 


By way of the iteration P m+1 (x) = (2J+3) (2J+1) " ' \ ' 5 ' jRmRm- 1 ■ ■ ■ RiRoPo(x) 
for h = 0 we obtain a generalised Rodrigues’ formula 


P m+ i(x) = (2m + 3)!! ^~ m/2+3/4 (^+4)- m/2+5/4 ^ ( a;3/2 ( a; + 4 ) 3/2 ^) m 1 ^ 3/4 (4 +x) 1 / 4 
where (2 n + 1)!! = (2 n + l)(2ra — 1) • • • 3 • 1. 

THEOREM 10. Let C n (x ) 6 e f/ie minimal polynomial of 2 cos and 0 n (x) the 
minimal polynomial of 2 cos (—). T/ien we have 



Pm(x~ 2) 

f=l 

if 2 m + 1 is a prime number, 


©2m+l (^) 

Z[x] 

otherwise. 




(4.26) 


Q m (x - 2) 

'=1 

= X 

Z[x] 

if m is a power of 2, 

if m is an odd prime number, 

otherwise. 




(4.27) 

Q2m-\-l{ % 2) 

(-l) m P m (-x 2 ) 

f=l 

if 2 m + 1 is a prime number, 

\m-\-2 (^) 

C4m+2 (x) 

€ Z[x] 

otherwise. 




(4.28) 

Qm( X - 2) 

(-l) m Q m (-x 2 ) 

f=l 

if m is a power of 2, 

1 ’ (4.29) 

otherwise. 

Cirri, (x) 

C^m (%) 

^ € Z[x] 

(-1 ) m P m (-x- 

2) _ I4i(x) _ 

f=l 

if 2 m + 1 is a prime number, 

C’ 2 m+l(x) 

C 2m+ i (x) 

€ Z[x] 

otherwise. 


COROLLARY. Let p n (x) be the minimal polynomial of 2 cos (^) — 2, r n (x) be 
the minimal polynomial of 2cos (^) — 2, and ip n (x) be the minimal polynomial of 
—2 cos — 2. T/ien 

Pm(x) = Pd 0*0, <2m(z) = II T 2 d(x), Pm{x) = <Pd(x), 

d\2m-\-l d\m d\2m-\-l 

d> 3 m/d is odd d> 3 


(4.31) 
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and 

Qm(x') — Pmi{x')'P mi (^x) — j | Pd{x)Td(xf), 

d\2mi~\-l 
d> 3 

when m = 2m\ is even. If m = 2m\ + 1 = 2 r {m r + 1) — 1 is odd, with m r 
even, so that for 1 < j < r — 1, rrij = 2mj + \ + 1, we have 

mi + 1 _ m 2 + 1 _ _ m r -\ + 1 _ 

777,2 + 1 W-3 + 1 777 r + 1 


then 


i —1 r —1 

Qm(x') — Qirir(x') || Qm.j+ 1(®) — Pm r .|-i (•®)^m T .+i (®) || Qmj +1 (^)> 
7=1 7=1 


(4.32) 


2?T7j,_|_7 


(4.33) 


(4.34) 


Proof of Theorem 10. We recall that the minimal polynomial of an algebraic number 
/3, is defined to be the rnonic polynomial of minimal degree, with rational coefficients, 
which has f3 as one of its roots. Such polynomials often exhibit structural properties, 
such as 4> n (x), the minimal polynomial of a primitive nth root of unity, e(k/n), with 
(k,n) = 1, which satisfies 

*»-i=n *-(*>■ (435) 

d\n 

It was shown in m that when n is a prime number p, then the minimal polynomial 
0 n (x) of 2 cos (2ir/n) is given by 0 n (x) = / n (x), with 


fn(x) = J^(-l)* 


k =0 


n — k 
k 


[^1 


n—2k 


E (- 1 ) 1 


k =0 


n — k 
k- 1 


n—2k-\-l 


and that 0 n (x) | f n (x) for all n € N. By algebraic manipulation we have P m (x — 2) = 
f m (x), and hence 114.261) . In fact, for p a prime number, we can write 0 p (2x) = 
2(p~i)/ 2 \I/ p (x), where v I' n (x') denotes the minimal polynomial of the algebraic number 
/3(n) = cos (2ir/n). 

It was shown by Watkins and Zeitlin (22] that analogous formulae to (|4.35l) for 
4' n (x) are given by 


T ni+1 (x) - T ni (x) = 2 ni || 4 r d (x), ?7 = 2m + 1 is odd, 

d\n 

T ni+ i(x) - T ni _i(x) = 2 ni ^ d (x), ?7 = 2m is even, 

d\n 

from which we can establish the explicit formula 


'J'n(x) 



(4.36) 

(4.37) 


so that deg 4 , n (x) = 1 if 77 = 1, 2 and <f{n)/2 if n > 3. From this one deduces that 
C'n(x), the minimal polynomial of 2cos7r/?7, is given by 


Ci = 24 / 2 


C n (x) = 2^'H 2r 


77 > 2. 
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It follows that deg C n (x) = 1 if n = 1 and 0(2n)/2 if n > 2, the zeros of C n (x),n > 2, 
are 2cos(7r k/n), with k = 1 1 and (k,2n) = 1. Hence each of the expressions 

in (14.271) . (14.281) . (14.291) and (14.301) are in Z[x] and equal to 1 for the prime conditions 
stated. In fact, for n an odd integer, we have Q n (—x) = (—1 )^ 2n P 2 C n (x). 

When n is a power of 2 we use the identity 2 2 ™ 2 'f> 2 rn (x) = 2T 2 m-2(x), and the 
result follows. 

For completeness we state the case that n is an odd prime power p m , with 
p = 2q + 1, for which we have 

2 p "‘- 1 ( p —D/2^ pm( x ) = 2 + L 

The Corollary can be deduced either from the relations (14.261) . (|4.27|> . (|4.28|) . (|4.29l) 
and (I4.30p . with x — 2 replaced by x, in conjunction with (13.181) . (13.191) and (13.201) . 
or directly from the properties of Fibonacci and Lucas polynomials. In particular, 
for p a prime number, the pth Fibonacci and Lucas polynomials are irreducible 
and so their roots are respectively 2 i times the real and imaginary parts of the pt\i 
cyclotomic polynomial, except for the root 0 in the Lucas polynomial case (see for 
example Koshy [8] 2001, p462). □ 

Example (of polynomial factorisation and special values). Taking m = 223 = 
2 5 (6 + 1) — 1, we have 

223 , ]c\ 

Q 223 (*e) = ^ ( 2 > I 1 ) x ' fc = p 3 {x)'P 3 {x)Q 7 (x)Qu{x)Q 28 (x)Q 5 ( i (x)Qn 2 (x). 

k =0 ' ' 

Without proof, some special values of the polynomials P m (x) and Q m (x) are 
Qm(-4) = 2(-l) m , P m (- 4) = (-l) m , Q m ( 0) = 2, P m (0)=2m + 1, 

Qm( 1) = £ 2771 ( 1 )) Pm{ 1) = £2m+l, Qm(l) = ^2m+2(l), Pm{ 1) = p 2m+l- 

5 On Minor Recurrence Relations 

The matrices M 0 (m,r) and M e (m,r), defined after the proof to Theorem 1, respec¬ 
tively generate our sequences and for 1 < j < m, via the recurrence 

matrix R m . In consequence, the i x i minors of these matrices sequence form a set 
of sequences in their own right, which also obey (different) recurrence relations, for 
2 < i < nr — 1. The 2x2 minors then correspond to the difference between consecu¬ 
tive convergents /^ m,k \ after multiplying through by the 

product of the two sequence terms which form the denominators of the convergents. 

We now briefly outline the general theories underpinning these ixi minor recur¬ 
rence sequence properties. For 1 < j < m. let the m sequences {yjk} < kLi be defined 
by an m x m initial value matrix, and an m-th order rational linear recurrence ma¬ 
trix K rn , of the form (13.291) with ao = 1, described in the proof of Lemma 3.1. Then 
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Ujk obeys the recurrence relation 

m 

Vjk = — (®i Dj (k— l) T a 2 Vj (k— 2 ) + • • • + CL m yj ( k—m )) > where K m {x) = ^ ^ a m —jX\ 

3=0 

with K m (x) the characteristic polynomial of K m . 

In this general setting we also assume that the system of polynomials K m (x) are 
orthogonal, and so satisfy a three-term recurrence, whose measure is supported on 
some interval [a, 6] € M. This ensures that the roots X m i, ■ ■ ■, X mm , of the polynomial 
equations K m (x) = 0, are distinct, real algebraic numbers lying in the interval [a, b] 
and that these roots interlace. Hence for m > n, there is a root of K m (x) = 0 
between any two roots of K n (x) = 0. 

We note that the condition ao = 1 produces a system of normalised roots, so 
that X m i x ... x X m m = 1. We also note that the minimal polynomials for each of 
the algebraic numbers X mj divides the characteristic polynomial K m (x). 

As stated in the proof of Lemma 3.1, the sequences = (1, X m j, A, 2 nj -, Af nj -,...), 
with j = 1,2,... ,m form a basis for the solution space for all possible sequences 
satisfying this recurrence relation, and so for any possible starting values. 

Let { 2/1 tefclfclo’ ■ ■ ■ ■ 0 * sequences generated by an initial value 

matrix and the linear recurrence, so that in matrix form we can write 

/ 2/10 yn 2/i 2 2/13 ••• \ 

2/2 0 2/21 2/2 2 2/2 3 • • • 

Y= . . . . . , 

\ ViO llil Vi 2 Vi 3 • • • / 

and consider the sequence formed by successive i x i determinants 


De 


Vie Dii+i yie+2 ■■■ 2/1 i+i—i 
2/2t 2/2t+i 2/2 1+2 ■■■ 2/2M-J-1 


£ = 0 , 1 , 2 , 3 ,... 


ya yu+i yu+2 ■■■ yu+i-i 


The space of all such sequences D{Y) = (Z?o, D 1 , D- 2 ,...) is spanned by the D’s that 
you get by choosing an array Y of the form 

/l 7i 7 1 7i ••• \ 

1 72 72 72 • • • 

^ . . . . 5 

V 1 7* li 7f •••/ 


where the 7 r are i distinct values chosen from the m distinct, real eigenvalues X m j, of 
the recurrence matrix. Here the order of the rows of Y is irrelevant as we are simply 
looking at how to choose an i element subset of an rn element set. It follows that there 
are (™) ways to choose such a Y , and from determinant theory the resulting D(Y ) is 
itself a geometric progression of the form (C, C A, CA 2 ,...), where A = 7172 • • • 7*. 
Hence, generically, the space of all such “z x i minor sequences” must be the solution 
space of a linear constant coefficient recurrence of order at most (™). We have just 
proved the following lemma. 



43 


LEMMA 5.1. For i > 1, the sequence formed by successive i x i determinants D^, 
as defined above, obeys a linear constant coefficient recurrence of order at most (” l ). 
If the eigenvalues of the minor recurrence matrix all have absolute value less than 
1, then the sequence of i x i determinants Di will converge to some number a. 


Example (or minor recurrence relation coefficients). We illustrate this lemma with 
the i minor recurrences corresponding to the denominator generating function Prj(x). 


{—5, -7,-4, — lj — nl 

f7 IQ 292 233 1223 356 63 T. 

V ’ ±y ’ 11 ’ 11 ’ 121 ’ 121’ 121’ 

(a _72 _63 _356 1223 233 292 

1 + 11’ 11’ 121’ 1331 \ 1J31 ’ 14641 


2 4 

1331 ’ 1331’ 
19 


1 


14641 


1 


r-, _± _ l __ J_ 

l 1 ’ 11’ 121’ 1331’ 

{~Tl } 


14641 


14641 ’ 


161051 ’ 1771561 


Here the recurrence coefficients for the sequences are given in the topmost 

entry, and those for the 2x2 minors the second from top entry, so that they satisfy 
a recurrence relation with recurrence polynomial 


+0 


—7x 9 +19x 8 —^^-x‘ + 

11 


233 1223 = 356 , 63 , 72 

- X° -arH- x 4 -arH-. 

11 121 121 121 1331 


■x 2 — - 


-x+ 


1 


1331 14641 


= 0. 


The sequences of m — 1 x m — 1 minors appear to have the recurrence polynomial 


yV (-ir 


1 /2m - k\ k 

2 m + 1 — 2 k\ k J 


so that when rn is odd the recurrence polynomial can be factorised as 


m 


n 



2 — 2 cos 


/ 27 vk \ 

^ 2771+1 J 


2 m + 1 


= 0. 


In general, the (™j + 1 term recurrence relation that our sequences of 2 x 2 minors 
obey can be used iteratively to obtain the coefficients Cjk(v ), such that 


/("y'l 

’ yr r+1 
As(m,k) 

t/ V+1 




l 


®-( m,k ) aj.( m,k ) 
J'r ^ r +l 


^r+1 
^r +1 


< 



2 ’ 


which is of a similar form to that of Dirichlet’s Theorem for standard continued 
fraction convergents. 
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Appendix 


Table of m-dimensional integer interlacing Fibonacci sequences for m = 

1, 2, 3,4, 5, with 1 < j < m, and 1 < r < 10. 


m 

j/r 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

i 

1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

2 

1 

2 

-1 

3 

-2 

7 

-5 

18 

-13 

47 

-34 


2 

1 

-1 

4 

-3 

11 

-8 

29 

-21 

76 

-55 

3 

1 

3 

-2 

2 

-17 

22 

-29 

269 

-357 

474 

-4406 


2 

2 

-3 

4 

-37 

49 

-65 

604 

-802 

1065 

-9900 


3 

1 

-2 

3 

-29 

39 

-52 

484 

-643 

854 

-7939 

4 

1 

4 

BO 

46 

-271 

1702 

-10855 

69499 

-445420 

2855494 

-18307378 


2 

3 


108 

-675 

4293 

-27459 

175932 

-1127763 

7230222 

-46355652 


3 

2 


116 

-755 

4859 

-31184 

199988 

-1282310 

8221661 

-52713260 


4 

1 

E9 

73 

-487 

3160 

-20332 

130492 

-836893 

5366170 

-34405885 

5 

1 

5 

-5 

11 

-32 

99 

-3415 

10744 

-33830 

106545 

-335575 


2 

4 

-10 

28 

-85 

265 

-9156 

28817 

-90746 

285805 

-900180 


3 

3 

-11 

35 

-110 

346 

-11982 

37734 

-118845 

374319 

-1178980 


4 

2 

-9 

31 

-100 

317 

-11002 

34669 

-109210 

343988 

-1083461 


5 

1 

-5 

18 

-59 

188 

-6535 

20602 

-64906 

204447 

-643954 


Table of negative index m-dimensional integer interlacing Fibonacci sequences 
corresponding to the renumbered Fleck number quotients obtained from the se¬ 
quence terms = n$(—2r + 1, (—r + j) (mod n)). Values for m = 1, 2,3,4, 5, 

with 1 < j < m, and —8 < r < 1, are given. 







































































































































